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1.0 SUMMARY 


This document describes the theoretical formulation and use of L218 (LOADS), a digital 
computer program to calculate load coefficient matrices. The program calculates load coef- 
ficient matrices utilizing the force summation technique. The particular field of application 
of the program is the calculation of load matrices for the analysis of airplane flight loads due 
to air turbulence (both random and discrete). 

The structural mode shapes, inertia data, and aerodynamic forces of the vehicle must be formu- 
lated outside L218 (LOADS) and read as program input data from either cards or magnetic 
files. The load coefficient matrices are then calculated and stored on magnetic files for use 
in other programs. The type of load coefficient matrices calculated are: 

• Translational and rotational accelerations, velocities, and displacements (AVD) 

• Aerodynamic panel loads on lifting surfaces (PLDS) 

• Net panel loads on lifting surfaces (NPLDS) 

• Shears and moments (VBMT) 

The sensor equations, used in the equations of motion when using active controls, have the 
same form as the AVD loads. Thus, AVD loads can be used to add sensor equations to the 
equations of motion. 


2.0 INTRODUCTION 


The computer program L218 (LOADS) was developed for use either as a standalone program 
or as a module of a program system called DYLOFLEX developed for NASA under contract 
NAS 1-139 18 (ref. 1). The LOADS program (L218) was designed to meet the DYLOFLEX 
contract requirements as defined in reference 2. These requirements specify the need for a 
program to calculate dynamic load coefficient matrices for use in calculating design loads 
and/or for use as sensors in active control analyses. 

The objective of this volume is to aid those persons wishing to use this program to calculate 
load coefficient matrices. To meet this objective, the following items are presented: 

• The theoretical formulation of the problem 

• The design and structure of the program 

• Instructions to use the program 

A sample problem is also included in this volume to aid the user in execution of the program. 
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3.0 SYMBOLS AND ABBREVIATIONS 


BS 

BBL 

WL 

c 

F 

FpL (q, q) 

FpL (tig) 

fg 

f^PLDS 

^NPLDS 

I 

[Jl] 

[J 2 I 

L 

m 

M 

M j, M2> 

M4,M5 


T 

q, q, q 
V 
Vt 
“ g 

SI 


Body station, length 
Body buttock line, length 
Water line, length 

Load coefficient matrix for the excitation function 
Force 

Response aerodynamic force 
Gust aerodynamic force 

Streamwise distance from the reference point first encountering a gust to 
the points encountering the gust later, length 

Aerodynamic panel loads 

Net panel loads 

Mass moment of inertia 

Mass and static moment matrix 

Static moment and moment of inertia matrix 

Load 

Lumped mass 
Moment 

Load coefficient matrices (nonaerodynamic) of the generalized coordinate 
displacement, rate, and acceleration, respectively 

Load coefficient matrices (aerodynamic) of the generalized coordinate rate 
and acceleration convoluted with the Wagner function 

Euler transformation matrix 

Generalized coordinate displacement, velocity and acceleration, respectively 
Shear 

True velocity 
Gust angle 

Klissner indicial lift growth function 
Wagner indicial lift growth function 
Spatial frequency (co/Vj) 

Load coefficient matrix with gust penetration for the excitation function 
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Modal displacement 

Modal slope in streamwise direction 

Modal slope 

Mode shape (includes 0 and <Pq). 
Indicates convolution 



4.0 ENGINEERING AND MATHEMATICAL DESCRIPTION 


4.1 AXIS SYSTEM 

All mode shapes used as input to the LOADS program are assumed to be in the local motion 
axis system of each node and the coordinate locations of the nodes for which modes are 
defined are in the reference axis system. The local axis system is related to the reference 
axis system by the rotation of the reference axis system about its three axes into the local 
axis system using the Eular transformation triad. The following illustrates these angles in 
their positive sense: 



That is, (+) dihedral, 0, is (+) rotation; (+) incidence, 0, is (+) rotation; (+) sweep, S2, is 
(-) rotation. 

For a more detailed description of the axis systems used in DYLOFLEX, the user is referred 
to appendix B of the DYLOFLEX summary document (ref. 1). 

4.2 SIGN CONVENTION 

The LOADS program requires the input data be consistent with the following sign conven- 
tions: 

Position (reference axis system) 
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Mode shapes (motion axis system) 



The force axis system is defined as a left-handed axis system for convenience in calculating 
inertia forces and moments and in calculating the work done by external forces. 

4.3 LOAD EQUATIONS 

The load equations are formed using the force summation method described in reference 3. 
The general form os the load equations is similar to that of the equations of motion presented 
in reference 4. The general form of the load equations without gust penetration is: 


lL| = [Mj] jq} + [M2] {q} + [M3] |q} + [M4] \q\ * <I> + [M5] jq} * «I> 

+ |C3lig.+ 

where: 


[Mj] [M2] and(M3] 
[M4] and [M5] 



Load coefficients matrices (usually nonaerodynamic)t for the 
generalized coordinate displacement, rate, and acceleration, 
respectively; 

Load coefficients matrices (usually aerodynamic) for the generalized 
coordinate rate and acceleration convoluted with the Wagner 
function 

Load coefficient matrix for the excitation function (usually gust 
angle) convoluted with the Kiissner function 




Wagner indicial lift growth function 


tOne exception to this is the load coefficient matrices due to the noncirculatory aerodynamic 
terms of strip theory aerodynamics. 
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KUssner indicial lift growth function 


q = Generalized coordinate displacement, velocity, and acceleration 

responses, respectively. 

= Gust angle 

Relating these matrices in a physical sense^_[Mj] , [M^] , and [M3] are associated with the 
load resulting from structural response. [M4] and [M5] are associated with the load resulting 
from aerodynamic response. 

With gust penetration, the excitation function of equation ( 1 ) is frequency dependent and is 
defined as: 

1^1 [0]^cos jfgj j-i [0] sin jfg}^ (2) 

where: 

= cj/Vj, spatial frequency 

I fg I = streamwise distance from point or points first encountering gust to the points 
encountering the gust later 

[$] = lifting panels contribution to gust loads at designated gradual penetration load 

stations 

The load equations program, LOADS, (L218) generates the[Mj] to [M5] matrices and the 
[^] matrix. Gust penetration is always considered to be present in DYLOFLEX, since the 
no-gust penetratmn problem is a special case of the gust penetration problem with only one 
gust panel, i.e., ['^] degenerates into |^[ and fg is a single scaler. 

Four different types of loads are available in LOADS, namely: 

• Transational and rotational acceleration, velocities, and displacement 

• Aerodynamic panel loads 

• Net panel loads 

• Shears and moments 

4.3.1 ACCELERATION, VELOCITY AND DISPLACEMENT (AVD) 

Acceleration, velocity and displacement equations can be used as either load equations or as 
sensor equations. The AVD equations contain both translational and rotational motions and 
are written as follows: 

Translational acceleration or accelerometer equation = SCALRl [0] jqj = [M3] (q| 

Translational velocity = SCALR2 [0] |q| = [M2] {q| (3) 
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Translational displacement 

= SCALR2 [<t>] jq [ 

= [Mil 


Rotational acceleration 

= SCALR3 [<Pg] jq[ 

= [M 3 ] 

|q‘l 

Rotational velocity or rate gyro equationw 

= SCALR3 [<f>g] Jqj 

= [Mo] 

|q}(3) 

Rotational displacement 

= SCALR3 [(l>g] jqj 

= [Mj] 



where: 

[<p] = Matrix of modal displacement 

[<j)g] = Matrix of modal slopes 

SCALRl, SCALR2, 

SCALR3 = Scaler multipliers 

The scaler multipliers are included only for user program flexibility, such as changing accelera- 
tion loads from inches/sec/sec to g’s or for changing rotational velocity from radians/sec to 
degrees/sec, etc. 


The AVD load coefficient matrices Mj, Mo, and M 3 are only a function of the mode shapes, 
with the direction of the loads being a function of the axis system used in defining the modes, 
All input mode shapes to DYLOFLEX are in the local motion axis system for each structural 
node. In order to obtain AVD loads in different axis systems from the local motion axis 
system, a transformation from the local motion axis system to some arbitrary axis system is 
performed by: 


and 







4>y 

= 

T 


<l>y 


arbitrary axis 





local axis 


■00 x' 




■ 00 x" 

00X 

= 

T 


00y 


arbitrary axis 



_ 00 Z_ 


axis 


(4) 


where [T] is the Euler transformation triad defined as: 


m = [X] [Y] [Z] 
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or some other combination of [X] , [Y] , [Z] and where: 


* 


1 

0 

0 

X 

= 

0 

COS0JJ 

sin0j^ 



0 

-sin0jj 

COS0JJ. 

- - 


COS0y 

0 

sin0y 

Y 

= 

0 

1 

0 

_ ^ 


-sin0y 

0 

COS0y 



COS02 

sin0^ 

0 

Z 

= 

-sin02 

COS02 

0 



0 

0 

1 


( 5 ) 


In addition, if the input mode shapes are not at the desired location for a load station 
(location where load coefficients are calculated), interpolation of the mode shapes to a load 
station is required. Interpolation on thin bodies can be performed using the interpolation 
array generated in INTERP (ref. 5). However, the interpolation methods used in INTER? 
are not directly applicable for calculating mode shapes at load stations on slender bodies. 
Consequently, two interpolation methods for use on slender bodies are available, one which 
utilizes a linear interpolation scheme between points and the other which is similar to the 
motion point interpolation method in INTERP. The latter method attaches rigid links from 
the node of known motion to the node where interpolation is required. 


To interpolate to a point x on a slender body, given two reference structural node points 
I and I+l and their coordinate locations (BS, BBL, WL): 


I X I+l 

A =^S^-BSjy(BSn.j -BSi) 

0Xx “ 

0Zx = A ^0ZI+1 - 0zij + 0ZI + LT ^00 Xx) 
*^0xx “ ^ (*^®xi+i ~^0xi^ "^^xi 

Wy, =A(*9yl+l-*«yl) 

00Zx. =A(^0zi+i-'^0zi) +<^0zi 

where: 

LT = BBL^-BBLj 
LTT= WLj^-WLj 
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To interpolate to a point x on a slender body, given only one reference structural node point 
I and its coordinate locations (BS, BBL, WL): 

I X 


•»Xx 


*Vx 

= 0yi-LB(002j) + LTT {00xx> 


= <Pzi + LB (<^>0yj) + LT (00xx ^ 

♦«Xx 

= <^0X1 

^9yx 

= «0yl 

Wxx 



where; 

LB = BS^ - BSj 

LT = BBL^ - BBLj 

LTT = WL^ - WLj 

The^interpolated modes are used to calculate the AVD load coefficient matrices Mj, M 2 , 
and M 3 in equation (3). 

4.3.2 AERODYNAMIC PANEL LOADS (PLDS) 


The aerodynamic panel loads are the loads produced on the panels of a thin body by the 
aerodynamics from the structural and gust response of the vehicle. The equation for the 
aerodynamic panel loads is: 


IFpLDsl = SCALR2 [FpL (q)]jq[ + SCALR2 [FpL (q) ]|ql * ^ 

+ SCALR3 (FpL(«g)] {og[ ♦ 'P 
or 

{FPLDS) = [M4] {q} + [M5] |q) + [?] jttg) 

where: 
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SCALR2,SCALR3 = Scaler multipliers 

FpL = Matrices of aerodynamic panel forces on thin bodies 


r = Wagner indicial lift growth function 
' = 1 if unsteady aerodynamics is used 
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( = Kiissner indicial lift growth function 
' = 1 if unsteady aerodynamic is used 

The coefficient matrices of aerodynamic panel forces used in the formulation of the aero- 
dynamic panel loads are generated in L2 1 7 (EOM) (ref. 4) in order to calculate the generalized 
aerodynamic forces. These aerodynamic panel force matrices are saved on file for use in 
L218 (LOADS). In this program, the user specifies the aerodynamic panels where the panel 
aerodynamic loads are to be calculated. This program will sort the aerodynamic p^el force 
matrices for these panels and assemble the appropriate matrices [M4] , [M5] and [^] . The 
scaler multipliers are included only for user flexibility in the program use, such as changing 
the resulting load units. 

4.3.3 NET PANEL LOADS (NPLDS) 

The modeling of a flight vehicle in DYLOFLEX for calculating dynamic loads requires the 
vehicle be represented as a number of structural and aerodynamic panels. In order to obtain 
net panel loads on thin bodies, it is necessary to sum aerodynamic and inertia forces on each 
panel. The expression can be simply stated as; 

jf^NPLDsj =SCALRj [m] [?J] |q} + jFpLosj 
or 

I^NPLDsj “ f^3^ {'dl + [M4] jq[ *4> + [M5] jq| *4>+[0] |o:g| (9) 


where: 


[m] 

= mass matrix 

SCALRl 

= scaler multiplier 


The most convenient panel representation of the vehicle would be that the structural and 
aerodynamic panels be identical. Practically, this usually is not the case, and consequently, 
interpolation of the aerodynamic panel forces may be necessary to match the structural 
panels. Two interpolation methods are included in L218 (LOADS) to interpolate the 
aerodynamic panel forces to the structural nodes (locations where inertia forces are calcu- 
lated). 

The first interpolation method requires a manual calculation of 1) which aerodynamic panels 
are to be used and 2) the percent of the aerodynamic panel forces to be added to each 
structural panel where the net panel force is required. This data is used to calculate a matrix 
called a weighting factor matrix, P, which is multiplied times the aerodynamic force matrices 
FpL (q), FpL (q), and Fp^ (oig). From this matrix calculation, it can be seen that the size 
of the weighting factor matrix must be the number of structural panels where net panel 
forces are required per surface by the total number of aerodynamic panel forces on this 
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surface. (The aerodynamic panel force coefficient matrices are generated in L217 (EOM) and 
include all of the aerodynamic panels on a surface.) As an example, to calculate the panel 
aerodynamic force at the structural node (which covers the area outlined by the dashed 
line) requires the aerodynamic panels 2, 3, 4, 6, 7, and 8 with only 0.5, 1.0, 0.25, 0.5, and 
0.5 of each aerodynamic panel forces, respectively (see fig. 1). 




Aerodynamic panel 
node 

Structural panel 
node 


Aerodynamic panel 
boundary 
Structural panel 
boundary 


Figure l.—Net Panel Load Example With Dissimilar 
Structural and Aerodynamic Panels 


A weighting factor matrix can be generated that would have dimensions 1 by 8 and would 
appear as; 


[P] = (0 0.5 1.0 1.0 0 0.25 0.5 0.5] 

The aerodynamic panel force applied to the structural panel would now be 

[P] [FpL(q,q,«g)] 

substituting this into equation (9) yields 

|FnPLDs1 =SCALR1 [M] [<t>] |q’l + SCALR2 (P] (FpL(q)] |q} * <I> 

+ SCALR2[P] (FpL(q)] {q} * 4> + SCALR3 [P] (FpLCttg)] jag[ * 'P 
or 

j^NPLOs) = t^3] jq| + [M4] jq) ♦ <P + [M5] |qj * <t> + [^] jbg[ * vp 

It is important to note that [P] in this example is a 1 by 8 matrix because the net panel load 
is only required at one structural node and the [m] would be only the mass for the structural 
node. If the net panel load is required at three structural nodes, the size of [P] would 
be a 3 by 8 matrix and the mass matrix would be a diagonal 3 by 3. In this program, and for 
this option, the user specifies the structural panels where net panel loads are to be calculated 
and the weighting factor matrix. The program will sort the structural panels to be used, 
calculate the inertia forces on these panels to form [M3] , calculate the aerodynamic forces 
using the weighting factor matrix to be applied to the aerodynamic at these structural nodes, 
and form [M4] , [Mj] , and (^] . The scaler multipliers, SCALRl, SCALR2, and SCALR3, 
are included only for user program flexibility, such as changing units. For example, if the 
mass matrix is actually input as weight, SCALRl could be used to change the resulting 
inertia force matrix M3 to the proper units. 
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The second interpolation method is similar to the previous method since it calculates new 
aerodynamic panel forces that are applied to the structural panels. This interpolation method 
requires as input the structural areas that each structural node encompasses and where net 
panel loads are required to be calculated. The method calculates the aerodynamic pressures 
applied to each aerodynamic panel node on a surface by dividing the aerodynamic panel force 
with the aerodynamic areas. These aerodynamic pressures at the aerodynamic nodes are 
interpolated to the structural nodes using the surface spline interpolation routine from L2 1 5 
(INTERP). These interpolated aerodynamic pressures are then assumed constant over the 
structural panel. Thus, to calculate the aerodynamic forces at the stnictural nodes, the 
interpolated aerodynamic pressures are multiplied by the structural areas. These new 
interpolated aerodynamic panel forces are then used in equation (9) for calculating net panel 
loads. It is apparent from this interpolation scheme that if the pressure gradient varies rapidly, 
the accuracy of the results degenerates rapidly. Thus, this interpolation method is only valid 
over regions where dp^/dx^ and dp^/dy2 are equal, or approximately equal, to zero. 

4.3.4 SHEARS AND MOMENTS (VBMT) 


The basic equations to calculate shears and moments at a point are, respectively: 
Shear = 2 forces 

Moments = 2 moment arms x forces + 2 moments 


( 11 ) 


The forces and moments of equation (11) consist of both structural inertia forces and moments 
and aerodynamic forces. 


The structural inertia forces at a structural node (i) where the mode shapes are calculated 
and in the force axis system oriented in the local axis system are given as: 


Fx 


m 

0 

0 

0 

-mz 

-my 

Fy 


0 

m 

0 

mz 

0 

-mx 

Fz 

i 

0 

0 

m 

my 

mx 

0 


00X 

00y 



( 12 ) 


= [Ji]i[0]i|q} 


where: 

m = Lumped mass at node i 

x,y,z = Distance from the structural node to the mass c.g. 
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The sign convention used for the arms (x,y,z) is: 


+ aft 

+ up • measured from the structural node to the mass eg location and in the 
local axis system 

+ outbd 


The moments at this structural node (i) and in the force axis system oriented in the local 
axis system are also given as; 





my I^x nixy 

-l 

Mx 


0 mz 

-mxz 

My 

= 

-mz 0 

mx mxy lyy 

myz 

[M.J 

• 

-my -mx 

0 -mxz myz 

Izz 

- 


•^x 

■ 

<t>ey 




( 13 ) 




The inertia matrices J j and J2 are the mass and rotary inertia data transferred from the 
inertia eg location to the structural node location. 


That is: 


XX 

^XXeg 

+ 

2 

my^ 

+ 

mz^ 


— 1 

+ 

2 

mx 

+ 

mz^ 

yy 

y^cg 





^zz 

II 

N 

N 

0 

CTQ 

+ 

2 

mx'^ 

+ 

my^ 


These structural inertia forces and moments in the force axis system are transformed from 
the local oriented axis to the inertia oriented axis system by an Euler transformation matrix. 


and 




( 14 ) 


( 15 ) 
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Where [Tjjjgj-^ja] is the Eular transformation matrix for the structural inertia forces and 
moments at structural node (i). 


The aerodynamic forces at aerodynamic node (i) in the force axis system transformed from 
the local oriented axis to the inertia oriented axis system are: 



similarly for Fpj^ (q) and Fp^ (^g) 

where [T^gj-Q^ is the Eular transformation matrix for the aerodynamic forces at aerodynamic 
node (i). 


To obtain shears and moments at any point, equations (12), (13), (14), and (15) are sub- 
stituted into ( 11 ) for the inertia loads and equation (16) is substituted into ( 11 ) for the 
aerodynamic loads and summed over the required nodes for the load in the force axis 
system oriented in the inertia axis system. The resulting equations for shear and moment 
are: 


|Shear}fQj.gg 2: ([Tinertial n,! [^.] {iif) 

+ 2 : [Taero! ([FpL ®) P) • 

+ [PpL W>1 1^1 • » + (FpL <“g) l{“g} • '•' ) O’) 

= [M 3 ] {q[ + (M4] |q}* 4 > + [M5] |q| :,<*> 


Moment force axis = 2 ([Moment arms] [T^ngrtia] [J)] [<^] 

+ [Tinertial 021 ® ) (ill 
+ 2 [Moment arms] [T^gj-g] 

([PpL (4)] {q| * + [FpL (q)] |q[ * (18) 

+ [FpL («g)] {oigl* 4^) 

= [M 3 ] |q[ + [M 4 ] |q[ * 4> + [M 5 ] |q| * 4> + [?]|ag|* 4^ 
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To calculate loads on the various airplane components, the most efficient technique is to sum 
the shear and moment loads on each component up to the intersection of its supporting 
component (defined as a dummy node), and then transfer the dummy load on to the second 
component and add it to the loads carried directly on that component. For example, the 
shear forces and bending moments at the root of the fin are applied as discrete forces and 
moments at the fin-body junction, and then summed along with the forces and moments 
carried directly on the body in calculating body loads (fig. 2). 



Figure 2.— Dummy Node Representation 

The VBMT module of the LOADS program is designed and coded to calculate shears and 
moments by modeling the vehicle as a number of distinct surfaces (as shown in fig. 2) whose 
loads are transferred from one surface to another by the use of dummy loads located at 
dummy nodes. The choice of the number and location of the dummy nodes is the user’s 
responsibility. 

The general form of the shear and moment equation is shown in equations (17)and (18). 
However, for user flexibility, equations (14) and (15), and (16) are multiplied by a scalar 
multiplier SCATS and SCALA, respectively, for each structural and aerodynamic node (i). 
The scaler multiplier is used in calculating the amount of structural inertia load and aero- 
dynamic load at each node to be used in the shear and moment calculation. SCATS and 
SCALA are made equal to one by the program for structural and aerodynamic nodes, 
respectively, to be included and equal to zero for the nodes not to be included, when using 
load panel summation card 6.9.2 (sec. 6.3). By using load panel summation card 6.9.3, 
SCATS and SCALA can equal any fractional percent for panels whose fractional loads are 
to be used in the shear and moment calculation. Card 6.9.3 requires that the values of 
SCATS and SCALA be manually calculated and input for each structural and aerodynamic 
nodes where fractional loads are to be included. 

In addition, for user flexibility, scaler multipliers SCALE 1, SCALE2, and SCALES can be 
used to multiply the load coefficient matrices M3, M4 and M5, and^, respectively, for each 
surface using card 6.4. These scalers can be used to change units. For example, if the inertia 
matrix is input in weight units, SCALE 1 can be used to multiply M3 by 1/g so that M3 will 
still be a force. 
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Referring to figure 2 to calculate loads at the load station, the shear and moment equations 
(17) and (18) respectively would appear as 

Shear = 2 aft fuselage forces + dummy load forces at the dummy node (19) 

Moment = 2 moments due to aft fuselage forces and moments + 2 moments 

due to dummy load forces and nloments at the dummy node (20) 

The summation of the forces and moments due to the dummy loads can also be multiplied 
by a scaler SCALED, which can equal any fractional percentage of the dummy load to be 
included in the calculation of the loads at a load station. SCALED must be manually cal- 
culated and input for each dummy load to be included in a load calculation using card 6.9.4. 

Of particular significance are two values of SCALED: -2, and 2. In load calculations, cases 
will arise, for example, when the loads from the horizontal stabilizer for both the left and 
right sides need to be included in calculating aft body loads. In a symmetric analysis, the 
force Fy and moments and on the left and right horizontal stabilizer will counteract 
each other, resulting in zero load, and the F^, and My will add to each other. Conse- 
quently, the program is designed so that if SCALED = 2 for a dummy load, Fy, M^^, and 
M^ are set equal to zero and Fj^^, F^, and My are doubled. Thus, when calculating symmetric 
aft body loads that include the horizontal stabilizer, only the right horizontal stabilizer is 
used as a dummy load and SCALED = 2 for this dummy load. Similarly, for an antisymmetric 
analysis, SCALED = -2 and F^^, F^, and My are set equal to zero and Fy, M^^, and M^ are doubled 
for the specific dummy load. 

A problem also exists when using aerodynamic slender bodies for loads calculations. The 
problem is that the aerodynamic modeling of slender bodies does not match the structure 
model (see fig. 3). The aerodynamic model is a straiglit tube, whereas the structure model 
will differ from this (i.e., in an upswept aft fuselage). The aerodynamic force locations are 
on the centerline of the aerodynamic model. These locations may not even exist on the 
structural model. The purpose of calculating shear and moments at a specific load station 
is for use in determining whether the structure is strong enough to carry the load. Conse- 
quently, the load stations are located on the structure and the external forces must be what 
the structure feels. It is obvious that if the aerodynamic forces are located at points not on 
the structure, the structure will not feel the proper loads. To adjust for the aerodynamic 
modeling, it is recommended that the aerodynamic force be moved vertically without pro- 
ducing a moment, to the centerline or elastic axis location of the fuselage or store (if the 
slender body is a store). This shifting of the z coordinates of the aerodynamic force locations 
on slender bodies can be performed by using card 6.6. 

Once the shears and moments have been calculated in the force axis system and oriented 
in the inertial axis system, the shears and moments can be transformed to any orientation 
using the inverse of the Eular transformation and the angles as input by card 6.7. The shears 
and moments in any orientation is given as: 
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Figure 3.— Vertical Correction for Fuselage 
Load Station 
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5.0 PROGRAM STRUCTURE AND DESCRIPTION 


The program is structured as a system of five overlays (fig. 4): 


Main Overlay 

(L2 18,0,0) 

Program L218vc 

Primary Overlay 

(L218,l,0) 

Program RGEN 

Primary Overlay 

(L2 18,2,0) 

Program AVD 

Primary Overlay 

(L2 18,3,0) 

Program NPLDS/PLDS 

Primary Overlay 

(L2 18,4,0) 

Program VBMT 



Read general 
input 


Calculate 

AVD 

(acceleration, 

velocity, 

and displacement) 


Calculate 
NPLDS (net 
panel loads) 
and PLDS 
(panel loads) 


Calculate 
(shears, moments, 
and torsions) 


Figure 4— Overlay Structure of L218 (LOADS) 

The (0,0) main overlay, L218vc, controls reading of general data cards by the primary over- 
lay (1,0) and calls the proper primary overlay to perform the execution requested. The main 
overlay does not read input cards. 

The (1,0) primary overlay, RGEN, reads the general input data and the module cards ($AVD, 
SNPLDS, SPLDS, or $VBMT) which select the primary overlay for execution. 

The (2,0) primary overlay, AVD, reads the AVD card and tape input data, performs the 
AVD calculations, and writes the AVDTAP input tape. 

The (3,0) primary overlay, NPLDS/PLDS, reads the NPLDS (or PLDS) card and tape input 
data, performs the NPLDS (or PLDS) calculations, and writes the NPTAP (or PTAP) out- 
put tape. 

The (4,0) primary overlay, VBMT, reads and VBMT card and tape input data, performs the 
VBMT calculations, and writes the LTAP output tape. 

The input/output files for the overlays are shown in figure 5. 
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Printed 
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report 
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Figure 5.— Input/Output of L218 Overlays 




















Although L218 (LOADS) serves as a module of the DYLOFLEX system, it can be operated 
as a standalone program. When the program is run by itself, it becomes the user’s respon- 
sibility to generate input data in the format required by L218 (see sec. 6.4 for file formats). 
When used as a module in DYLOFLEX, L2 18 will receive input from L215 (INTER?) and 
L217 (EOM) (refs. 5 and 4, respectively). 


5.1 PROGRAM RGEN 

The first overlay called into execution by L2I8 is always the 1,0 primary overlay named 
RGEN. Program RGEN, overlay (L218,l,0) reads and interprets general input data from 
cards and determine the type of execution to be performed. RGEN will also perform the 
following optional operations: 

• Rename any of the input or output files. 

• Read mass data from cards and write a diagonal mass matrix file or nondiagonal J 
matrix file for use in overlays (L2 18,3,0) and (L2 18,4,0). These data files may also 
be saved for subsequent use in L218 (LOADS). 

• Redefine the order of rotation of the Euler angles. 

5.2 PROGRAM AVD 

Program AVD (acceleration, velocity, and displacement), overlay (L21 8,2,0) processes 
modal deflections (0), geometry data (BS, BBL, WL, 6^, dy, ) at required node numbers, 
and scalers to generate the appropriate coefficient matrices for use by L219 (EQMOD) or 
L221 (TEV156) (refs. 6 and 7, respectively). The coefficient matrices are used to calculate 
loads consisting of translational and/or rotational accelerations, velocities, and displace- 
ments at selected points on the structure. 

The axis system which defines the direction of the AVD loads may be either in the local 
motion axis system (the axis system in which the input mode shapes are defined), the 
inertia axis system (the axis system defined relative to the local motion axis system by the 
angular data d^, dy, and 6^ in the geometry input from INTERP), or any arbitrary axis 
system (an axis system defined relative to the reference axis system by 0^, dy, and d^ 
from card input). The transformation of the modal data from the local motion axis to the 
inertia or arbitrary axis systems uses the Eular transformation triad and requires that each 
mode-shape (/>^, <py, <t>^, (pQ ,(pQ , and 0^ exists at each required node. If any one of these 
mode shapes do not exist i?i the^input dafe, they are assumed to have zero mode shape. 

The modal data can also be interpolated to arbitrary locations. When the interpolation 
option is selected, axis transformation (rotational) is not allowed. Interpolational is per- 
formed in the local axis system to any point specified by reference axis coordinates. 
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There are three interpolation methods available, one method for thin bodies, and two meth- 
ods for slender bodies. 

The interpolation method for thin bodies uses the interpolation method from INTER? for 
the specified thin bodies. The mode shapes (p^, (py, and <Pq are assumed zero. Mode shapes 
‘Pz’ *^0 > “^0 interpolated using the appropriate SA array from INTER? and the results 
multiplfed b^-1 to retain the proper sign convention of the modes. The first interpolation 
method for slender bodies uses a linear interpolation method between two given reference 
structural nodes or from one given structural node to any node using rigid links (see sec. 
4.3.1). This method is suitable for elastic axes which are straight and preferably parallel ■ 
to the reference axis. The second interpolation method for slender bodies uses a linear inter- 
polation method which interpolates to a node using rigid links attached to a reference node. 
The equations are the same as the first method; however, the length of the rigid links LB, 

LT, and LTT and the reference structural node number (I) is input from cards. 

5.3 PROGRAM NPLDS/PLDS 

Program NPLDS/PLDS, overlay (L2 18,3,0) reads the specific card input data for net panel 
loads or aerodynamic panel loads. NPLDS (net panel loads) uses the m^dal deflection 
with the mass matrix and scaler multiplier to obtain the inertia forces (M3 matrix) on each 
surface. It also uses the aerodynamic force matrices generated in L217 (EOM) (ref. 4) 
from which selected nodes are extracted and multiplied by a scaler to form the M4, M3, and 
^ matrices. Only the modal deflection is used since net panel loads (also aerodynamic 
panel loads) are loads that act perpendicular to a lifting surface (thin body). PLDS 
(aerodynamic panel loads) is identifical to NPLDS except that the [M3] matrix is omitted. 

A net panel loads run will result in loads coefficient matrices written on file (NPTAP) and 
a panel loads run will result in load coefficient matrices written on file (PTAP). 

If the aerodynamic and structural panels are not identical when calculating net panel loads, 
two interpolation methods are available that interpolate the aerodynamic panel forces to the 
structural nodes. One method requires a weighting function matrix be manually calculated 
and input while the second method requires input of the structural areas corresponding to 
the structural nodes. 

For each requested surface and the selected sfructural nodes, the program calculates [M3] , 
writes it on NPTAP, forms [M4I , [M5] and [^] for all frequencies (k), and writes them 
on NPTAP (also on PTAP if requested). For a more detailed description of the computer 
program steps in calculating [M3] , [M4] , [M3] , and [?>] , the user should refer to volume 
II of this document. 


5.4 PROGRAM VBMT 

Program VBMT, overlay (L2 18,4,0) reads specific card input data and calculates shears and 
moments coefficient matrices for use in subsequent programs. All forces and moments are' 
transformed to the inertia axis system before shears and bending moments are calculated. 
After all calculations are performed, the shears and moments are rotated to the desired 
orientation. 



mil mill 


The calculation of the M3 load coefficient matrix, which consists of shears and moments due 
to inertia forces and moments, requires the mode shapes from the SATAP_generated in L21 5 
(INTERP) and the J matrix input in LOADS. The calculation of the M4, M5, and ^load 
coefficient matrices requires the aerodynamic force coefficient matrices Fp^ (q), (q), a„ 
from the EOMLOD tape generated in L217 (EOM). 

I^r eachjoad set and each requested surface and selected load, the program calculates [M3] , 
[M4] , [M5] and i^] for all frequencies (k) and writes them on LTAP. 

For a_detailed description of the computer program steps in calculating [M3] , [M4] , [M^] , 
and [^] , the user should refer to volume II of this document. 
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6.0 COMPUTER PROGRAM USAGE 


The program was designed for use on the CDC 6600. The machine requirements to execute 

L2 18 (LOADS) are: 

Card reader — To read control cards and card input data. 

Printer — To print standard output information, optional intermediate results and 

diagnostic messages. 

Disk storage — All magnetic files not specifically defined as magnetic tapes are assumed 
to be disk files used for input, output, and temporary data storage. 

Tape drive — For “permanent” storage of data. Magnetic files are copied to and from 
magnetic tapes with control cards before and after program exeuction. 

The program L218 (LOADS) is written in FORTRAN and may be compiled with either the 

RUN and FTN compiler. L218 may be executed on either the KRONOS 2.1 or NOS oper- 
ating system. 


6.1 CONTROL CARDS 

The following list is a typical set of control cards used to execute L2 1 8 (LOADS) using the 
absolute binaries from the program’s master tape. 

Job Card 
Account Card 


REQUEST (MASTER, F=l, LB=KL, VSN=66XXXX) 
REWIND (MASTER) 

SKIPF (MASTER) 


Retrieve the 
program from its 
master tape 


L218 


Prepare optional 
input data files 


Execute L218 (LOADS) 
Save optional output 
data files. 
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EXIT. 

DMP (0, field length) 
REWIND (lOUT) 

COPYBF (lOUT, OUTPUT) 

— End-of-record 


Card Input 
Data 


— End-of-File 


The following list is a typical set of control cards used to execute L218 (LOADS) using the 
relocatable binaries from the program’s master tape. 

Job Card 
Account Card 


REQUEST (MASTER, F=I, LB=KL, VSN=66XXXX) 
REWIND (MASTER) 

SKIPF (MASTER, 2) 

COPYBF (MASTER , REL218) 

RETURN (MASTER) 


LOAD(REL218, DYLIB) 
NOGO. 

L218. 


Retrieve the 
program from 
its master tape 


Prepare optional 
input data files. 

Retrieve DYLOFLEX 
alternate subroutine 
Ubrary, DYLIB. 

Load and execute 
L2 18 (LOADS) 


EXIT. 

DMP (0, field length) 
REWIND (lOUT) 

COPYBF (lOUT, OUTPUT) 
—End-of-record 

Card Input 
Data 

— End-of-file 
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6.2 RESOURCE ESTIMATES 


Field Length 

The core required for L218 is dynamic and changes with each input case. The following 
guidelines are given for each module: 

AVD — Required core is approximately 1 1 6000 octal 

+ (9) (NODES) + (6) (NODES) (MXNOE) 

where: 

NODES is the number of nodes requested and MXMODE is the 
number of modes from the SATAP tape. For thin body inter- 
polation, an additional core is required which is the larger of 2 (NODES) 
+ (NODES+3) (MXMODE) 
or (NODES+3) (NODES+4) 

NPLDS/PLDS — Required core is approximately 130150 octal 

+ (NODES) (6+MXMODE) + (NAERO) (3+NLOAD+MXMOE) + NK 

where: 

NODES and MXMODE are as defined above, NAERO is the number 
of aero panels from the equations of motion tape, NLOAD is the 
number of structural loads, and NK is the number of frequencies. 

VBMT — Required core is approximately 122400 octal 

+ (60) (NK) + (5) (NAERO) + (NODES) ( 1 8+ (6) (MXMODE) ) 

+ (MXLOAD) ( 1 0+(9) (NK) + (2) (NAERO) ) 

where: 

NK, NAERO, NODES, MXMODE are as defined above, and MXLOAD 
is the number of loads for a given load-set. 

Each module computes the total core requirement for a given run and prints it in octal as: 
REQUIRED DYNAMIC STORAGE= 

Thus, for similar problems, core can be set equal to the maximum required after the first 
run. 

If core is insufficient, a probably indication wUl be a CPU-01 error (address out of range). 
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6.3,1 PROGRAM RGEN 


The primary overlay RGEN reads program directive cards to: 

1 . Assure that the data being read is intended for L2 1 8 (LOADS) 

2. Determine which section of code (primary overlay) of L218 is to be executed next 

3. Determine what data and results are to be printed 

4. Determine input and output magnetic file names 

Program RGEN reads and processes card sets 1.0, 2.0, 3.0 and 7.0. RGEN also reads the 
“$” card of card sets 4.0, 5.0 and 6.0 and sets the control flag (KMOD) to cause L218 to 
execute the appropriate overlay to process the rest of that card set. 

Card sets 1 .0, 2.0 and 7.0 are simply the cards SLOADS, STITLE and SQUIT. The order 
of input for card set 3.0 is displayed in figure 7. 



Figure 7.— Flow of Card Set 3.0 input Data 
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Card Set 1.0 — Introduce L218 (LOADS) Card Input Data 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$ LOADS 

AlO 

Keyword introducing LOADS input. 


Card Set 2.0 — Title Card (Optional) 

Repeat the title card as desired to label the printed output. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

I-IO 

STITLE 

AlO 

Keyword identifying title card. 

11-70 

Title 

6 AlO 

Job title which can provide description 
of the job. 


Card Set 3.0 — Introduce General Loads Input Data (Optional) 

. ■ ) 
I 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$GEN 

AlO 

Keyword introducing general loads input. 
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Card 3.1 — File Names (Optional) 

Input names of files (tape or disk units) to be used for input/output data. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

IFILE (1)* 

AlO 

Keyword indicating the file to be named. 

11-20 

INAME (1) 

AlO 

Name given to the file. (Maximum of 6 characters) 

• 

• 

• 

• 

• 



41-50 

IFILE (3) 

AlO 

Keyword indicating the file to be named. 

51-60 

INAME (3) 

AlO 

Name given to the file. 




Use additional cards if necessary. 


• NOTE: 


IFILE (Keyword) 


DESCRIPTION OF FILE 

AVDTAP 

AVDTAP 

Final AVD output tape. 

NPTAP 

NPTAP 

Final NPLDS output tape. 

PTAP 

PTAP 

Final PLDS output tape. 

LTAP 

LTAP 

Final VBMT output tape. 

EOMLOD 

EOMLOD 

Equations of motion input tape for LOADS. 

SATAP 

SATAP 

INTERP input tape for LOADS. 

MASSTP 

MASSTP 

Diagonal Mass matrix input tape for LOADS. 

JTAPE 

JTAPE 



J-Matrix input tape for LOADS. 


Card 3.2 — Introduce Mass Data (Optional) 

The mass data is input on card 3.2.2 for surfaces defined on card 3.2.1. 

The mass data may be saved on tape MASSTP (card 3.1) for subsequent usage. 



KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1 

O 

MASS 

AlO 

Keyword introducing the mass matrix data. 
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Repeat cards 3.2.1 and 3.2.2 for each surface until all surfaces have been input. The surfaces 
must be input sequentially beginning with surface 1 and continuing until all surfaces have been 
input. A null surface should be input with zero’s. 

Card 3.2.1 — Mass Data Surface Definition 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

Surface 

AlO 

Keyword to introduce a surface's mass matrix. 

11-15 

IS 

15 

Integer defining the surface number of the mass 
matrix which follows. 


Card 3.2.2 — Diagonal Mass Matrix Elements 

Repeat card 3.2.2 until all elements for surface IS is defined. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

XMASS (1) 

7E10.0 

Diagonal mass matrix elements for surface IS (as 

11-20 

XMASS (2) 


defined in Card 3.2.1) in weight or mass units. 

21-30 

XMASS (3) 



31-40 

XMASS (4) 



41-50 

XMASS (5) 


• 1 

51-60 

XMASS (6) 



61-70 

XMASS (7) 




Card 3.3 — Introduce Inertia Data (Optional) 

The inertia data is input on cards 3.3.3 and 3.3.4 for surfaces defined on card 3.3.1. The 
inertia data may be saved on tape JTAPE (card 3.1) for subsequent useage. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

JMAT 

AlO 

Keyword introducing the J-matrix data. 
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Repeat cards 3.3.1, 3.3.2, and either 3.3.3 or 3.3.4 for each surface until the inertia data for 
all surfaces are defined. The surfaces must be input sequentially beginning with surface 1 
and continuing until all surfaces have been input. If a surface is omitted, it is assumed to be 
null. 

Card 3.3.1 — Inertia Data Surface Definition 


COLS. 

KEYWO^/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

Surface 

AlO 

Keyword to introduce a surface’s J-matrix. 

11-15 

IS 

15 

Integer defining the surface number of the J-matrix 




which follows. 

K, !! 

, 

1 ' 

U> 

o 

NODes 

AlO 

Keyword to introduce the number of structural nodes. 

[31-35 

INODE 

15 

The number of structural nodes in this surface. 


Card 3.3.2 — Inertia Data Type and Format 

Repeat cards 3.3.2, and 3.3.3 or 3.3.4 as necessary to read all components of the J-matrix 
inertia data for surface IS (card 3.3.1). 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

i-10 

JTYPE 

AlO 

Keyword identifying the type of J-matrix data to follow. 
The keyword must be one of the following: MASs, MX, MY, 
MZ, IXX, lYY, IZZ , MXY, MXZ , MYZ . 

The order must be followed as listed. However, if 
one or more is null, it is omitted. 

11-20 

MATYP 

AlO 

Keyword FtILl or SPArse indicating a full diagonal matrix 
or a sparse diagonal matrix. This will determine the format 
for the matrix data which follows: Card 3.3.3 or 3.3.4. 


The J-matrix can be partitioned into a number of diagonal matrices (MASs, MX, MY, 
MZ, IXX, etc.). These diagonal matrices may be full or sparse. 
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Card 3.3.3 — Full Diagonal J-Matrix 

Omit this card if the diagonal J-matrix is sparse (MATYP=SPArse on card 3.3.2). 
Repeat this card until INODE elements have been read for surface IS (card 3.3.1 ). 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

XJMAT (1) 

7E10.0 

J-matrlx data for JTYPE matrix defined in Card 3.3.2. 
This is a full diagonal matrix format (designated by 
MATYP in Card 3.3.2). 

• 

• 

• 

61-70 

XJMAT (7) 


Card 3.3.4 — Sparse Diagonal J-Matrix 

Omit this card if the diagonal J-matrix is full (MATYP=FUL1 on card 3.3.2). 
Repeat this card as necessary to read all desired nodes on surface IS (card 3.3.1). 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

NODE (1) 

15 

Structural node number (I) 

11-20 

XJMAT (1) 

ElO.O 

J-matrix value for node (I) and JTYPE matrix defined 




in Card 3.3.2. 

• 

• 

• 

• 

• 

• 

• 

« 

• 

« 

• 

41-45 

NODE (3) 

15 

• 

51-60 

XJMAT (3) 

ElO.O 

• 
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Card 3.4 - Euler Angle Rotation Order (Optional) 
Transformation order for AVD and/or VMBT (see section 4.3). 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

TRANsform 

AlO 

Keyword introducing the rotation order of the Euler 
transformation matrix for AVD and/or VBMt. 

11-20 

lOROER 

AlO 

Keyword XY2, XZY, YXl, YZX, ZXY, or ZYX, indicating 
the rotation order of the Euler transformation matrix 
(See Section 4 ) . 

(Default « XYZ) 

XYZ is the order of the matrix multiplication which 
means rotation about Z first, then Y, and lastly X. 


Card 3.5 — Maximum Number of Surfaces (Optional) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

MAXSUR 

AlO 

Keyword introducing a non-standard value of ISMAX. 

11-15 

ISMAX 

15 

New value for ISMAX (maximum number of surfaces) . 
(Default = 20) 


ISMAX may be changed as required. A larger value of ISMAX will result in more dynamic 
storage being used. 
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Card 3.6 — Temporary File Buffer Size (Optional) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

LBUF 

AlO 

Keyword introducing & non-standard buffer size. 

11-15 

LBUF 

15 

New value for LBUF. This is the buffer size used 
by FETADD for temporary files. 

(Default » 100) 


For large problems requiring excessive disk storage, it may be desirable to increase LBUF. 
This would directly increase the dynamic storage required but may significantly decrease 
the number of required disk accesses. 


Card 3.7 — End of General Card Input 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SEND 

AlO 

Keyword indicates the end of Card Set 3.0. 


Program RGEN reads the following load option card and transfers control to the main overlay 
which calls the proper overlay. The following are the card sets for the corresponding load 
options. 


Card Set 

Keyword 

Load Option 

4.0 

SAVD 

Accelerations, velocity, and displacements 
(Requests overlay (L218, 2, 0) 

5.0 

SNPLDS 

Net panel loads/panel loads (Requests over- 
lay (L2 18, 3,0) 

5.0 

SPLDS 

Panel loads (Requests overlay (L218, 3, 0) 

6.0 

SVBMT 

Shears, bending moments, and torsion ( 
(Requests overlay (L218, 4, 0) 


Program RGEN also reads the SQUIT card (card set 7.0) and transfers control to the main 
overlay to terminate the program execution. 
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6.3.2 PROGRAM AVD 


Card set 4.0 contains instructions and data required to form load coefflcient matrices used 
to calculate accelerations, velocities, and displacements. These load coefflcient matrices 
may be used as sensor equations in L219 (EQMOD) (ref. 6), or as load equations in L219 
(EQMOD), and/or L221 (TEV156) (ref. 7). Each of the load coefflcient matrices cannot 
exceed a size of 100 x 70 (100 loads, 70 degrees of freedom). If more loads are required, 
the user must repeat card set 4.0. The $AVD card terminates the previous load set (if any) 
and begins another load set. The final output file has successive load sets separated by 
end-of-flle marks. 

The order of input for card set 4.0 is displayed in figure 8. 

Omit card set 4.0 if no acceleration, velocity, or displacement coefflcient matrices are 
required. 



Figure 8.— Flow of Card Set 4.0 input Data 
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Card Set 4.0 - AVD Load Option 

Request AVD execution and introduce the AVD card set. 

Repeat card set 4.0 for each load set when multiple load sets are required (limit of 100 loads 
per load set). 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SAVD 

AlO 

Keyword requesting the AVD path of the LOAD Module 
and introducing the AVD card input. 

11-20 

print 

or 

no-print 

AlO 

Visual keyword*indicating that a print key may appear 
in Col. 21-23, 31-38 or 41-48. If printing is desired a 
keyword must appear in the appropriate column (s). 

21-30 

HAP 

AlO 

Keyword requesting the general print summary. 

31-40 

MATrices 

AlO 

Keyword requesting the general print summary and the 
output matrices as placed on AVDTAP. 

41-50 

CHEckout 

AlO 

Keyword requesting the general print summary, the 
output matrices as placed on AVDTAP, and intermediate 
calculations. 


*Note: All visual keywords are optional and are not read by the program. 


Repeat cards 4.1 -4.8 as required. 
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If 


Card 4. 1 — Surface Definition 

The surface card (Card 4.1) may be repeated with the same surface number IS if necessary. 


COLS. 

KEYWORD/ 

VARIABLE 

= 1 

FORMAT 

DESCRIPTION 

1-10 

Surface 

AlO 

Keyword to introduce the surface number. 

11-15 

IS 

15 

Surface number requested (a default maximum number of 
20, see Card 3.5 for increase). 

16-20 


5X 

Not Used. 

21-30 

CARDS * 
or 
TAPE 

AlO 

Keyword indicating the source of the mode shapes. 
Keyword CARDS means the mode matrices will be read 
from cards (Card 4.8). 

Keyword TAPE means the mode matrices will be input 
from the file SATAP. 

Default ; TAPE 

31-35 

I ROWS 

15 

If Col. 21-30 contains CARDS, Irows is the number of 
rows in the mode matrix (number of nodes on surface IS) . 
IROWS is left blan)c if Col. 21-30 contains TAPE. 

36-40 

I COLS 

15 

If Col. 21-30 contains CARDS, ICOLS is the number of 
columns in the mode matrix (number of modes defined) . 
ICOLS is left blank if Col. 21-30 contains TAPE. 


*See note following card 4.4. 


Card 4.2 — Comment (Optional) 

Repeat card 4.2 as desired to label the card data and printed output. 


COLS. 

1-2 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

C_ 

A2 

Keyword indicating a comment card that will be printed. 

3-80 


A8,6A10 

Available for comments. 
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Card 4.3. — Angular/Translational Matrix Selection Card 
Deflnes the type of AVD load and M matrix to be generated. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

A-r 

6(A2,1X 

Keyword , 

where A is one of the following for 



A3,4X) 

surface IS 

11-20 

A-F 


TA - for 

translational acceleration matrix 




TV - for 

translational velocity matrix 


• 


TD - for 

translational displacement matrix 

• 

• 

• 

• 


RA - for 

angular acceleration matrix 




RV - for 

angular velocity matrix 

51-60 

A-F 


RD - for 

angular displacement matrix 




and where 

F is one of : 




XYZ, XY, 

XZ, YZ, X, Y, Z where 




X, Y or Z 

defines the X, Y or Z mode component to be used 




in generating the matrix. 




For example: TA-X would be a keyword to generate the 




matrix for translational acceleration in the x direction 




using the 

mode deflection. RV-XYZ would be a keyword 




to generate the matrix for angular velocity in the 




X, y, z directions using the mode shapes. . ^ . 
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Card 4.4 — Surface Options 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-iP 

NODE * 

AlO 

Keyword indicating which nodes of surface IS will be used 


or 


to calculate the type of load defined on Card 4.3. 


NODE-ALL 


Keyword NODE means only specific structural node numbers 


or 


will be read by Card 4.5. Keyword NODE-ALL means all 


COORD 


structural nodes defined in INTER? (L215) will be used 
(Omit Card 4.5). Keyword COORD means nodal coordinates 
(X,y,Z) will be read on Card 4.7. 

Default : NODE 

11-20 

LOCAL * 

AlO 

Keyword indicating the resultant load axis system (See - 


or 


Section 4.0). 


INERT 


Keyword LOCAL means the local axis option is selected. 


or 


Keyword INERT means the inertia axis option is selected. 


ANGLES 


Keyword ANGLES means arbitrary rotation angles 6y» 6^ 

i 



will be read with Card 4.6. 




Default : LOCAL 

21-30 

SCALRl 

ElO.O 

Scale factor on TA Default: SCALRl =1.0 

31-40 

SCALR2 

ElO.O 

Scale factor on TV and TD Default: SCALR2 =1.0 

4i-50 

SCALR3 

ElO.O 

Scale factor on RA, RV, RD Default: SCALR3 =1.0 

51-60 

PRINT 

AlO 

A keyword to define the sequential print option after 


or 


each surface calculation. 


NO-print 


Keyword PRINT prints surface number, (X, Y, z) 




reference coordinate, node number, type of load. 




and matrix coefficients. 




Keyword NO indicates no printing unless selected on 




$AVD (Card Set 4.0) 




Default : PRINT 


*See note next page. 




Note: The following table summarizes valid combinations of keywords from cards 4. 1 and 
4.4. If the user chooses an invalid combination of keywords, the user will get fatal 
error diagnostic. 


VALID KEYWORD COMBINATIONS 


Card 4 . 1 

Card 4 . 4 

CARDS 

COORD 

LOCAL 

TAPE 

NODE 

ANGLES 

LOCAL 

INERT 

NODE-ALL 

ANGLES 

LOCAL 

INERT 

COORD 

LOCAL 


Card 4.5 — Structural Node Numbers, (Optional) 

This card is included only if the first keyword on card 4.4 is NODE or blank. 
Repeat card 4.5 as necessary until all structural nodes to be used as AVD loads. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

I INODE (1) 

1415 

Numbers of nodes (monotonic-increasing) specifically 




requested to calculate the AVD loads requested on Card 4.3. 




For example: 1357 would include nodes 1# 3, 5, 7. 

• 

• 

« 

• 


If a string of nodes is required in addition to some 

• 

• 


selected nodes, it might be requested by 




13 5 -10 13 14 where 5 -10 is interpreted to mean 




nodes 5, 6, 7, 8, 9, and 10. 

66-70 

IINODE (14) 
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Card 4.6 — Transformation Angles (Optional) 


This card is included only if the second keyword on card 4.4 is ANGLES. 
Repeat card 4.6 for all structural nodes retained. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

THETAX 

3E10.0 

Thetax(0x^/ Thetay(0yj, ThetazCe^^) defines a set of 

11-20 

THETAY 


axis rotation angles for node i. The order of rotation is 

21-30 

THETAZ 


as specified in Card 3.4. 


Card 4.7 - Interpolation* (Optional) 

This card is included only if the first keyword on card 4.4 is COORD. 
Repeat card 4.7 for each coordinate location desired. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 


DESCRIPTION 

1-10 

X 

ElO.O 

BS 

Coordinates in the reference axis 

11-20 

Y 

ElO.O 

BBL 

• system of the point at which AVD 

21-30 

Z 

ElO.O 

WL 

loads are to be calculated. 

31-40 

ALB 

ElO.O 

LB (interpolation coefficient) = BS-BS.,_,,_ 

XbUB 

Omit if ISUB = blank 

41-50 

ALT 

ElO.O 

LT (interpolation coefficient) = BBL-BBL^e™ 

XoUo 

Omit if ISUB = blank 

51-60 

ALTT 

ElO.O 

LTT (interpolation coefficient) = WL-WL^c™ 
Omit if ISUB « blank 

61-65 

ACODE 

A5 

Interpolation option code*, 
blank or PAR 
Omit if ISUB > 0 

66-70 

ISUB 

15 

Reference structural node number (ISUB) 


*Interpolation Methods (see section 4.3,1). 
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ACODE/ISUB 

METHOD 

APPLICATION 

blank/blank 

Interpolation method used 
in INTERP tL215) 

Thin bodies 

PAR/blank 

Linear interpolation between 
structural nodes or extra- 
polation from a structural 
node. 

Slender bodies with elastic 
axis parallel to reference axis. 

blank/ref . node 
no. ISUB 

Interpolates to a node using 
rigid links of lengths ALB, 
ALT, and ALTT attached to the 
reference structural node ISUB 

Slender bodies 


Note; For cases using LB, LT, LTT, a particular reference node (ISUB card 4.7) can be 
used only once for a surface calculation. If more than one reference to the same 
node is required, the surface data must be repeated for each additional reference 
(card 4. 1 ). 


Card 4.8 — Card Input of Mode Matrix (Optional) 

This card is included only if the keyword on card 4. 1 is CARDS. 

Card 4.8 is repeated only for each mode matrix required, as determined by card 4.3, and 
in the order [<l)^] , [0y] , [0^] , [4>d-^] , l<t>dy] , [<Pdz] ■ 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

(1) 

7E10.0 

i^^ row of the mode shape matrix. 




The mode shape matrix will be read by rows, repeating 

61-70 

• 

4> (7) 


this card as necessary. Each row begins as a new card. 
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Card 4.9 - End of AVD Card Input 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SEND 

AlO 

Keyword indicates the end of Card Set 4.0. 
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6.3v3 PROGRAM NPLDS/PLDS 


Card set 5.0 contains instructions and data required to form net panel loads or aerodynamic 
panel loads using data read from the Interpolation (L215) tape (SATAP) (ref. 3), and the 
equations of motion (L217 tape EOMLOD) (ref. 5). 

Each of the load coefficient. matrices cannot exceed a size of 100 x 70 (100 loads, 70 
degrees of freedom). If more loads are required, the user must repeat card set 5.0. The 
SPLDS or SNPLS card terminates the previous load set (if any) and begins another load 
set. The final output file has successive load sets separated by end-of-file marks. 

The order of input for card set 5.0 is displayed in figure 9. 

Omit card set 5.0 if no net panel load or panel load coefficient matrices are required. 



Figure 9.— Flow of Card Set 5.0 Input Data 
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Card Set 5.0 — NPLDS/PLDS Load Option 


Repeat card set 5.0 for each load set when multiple load sets are required (limitation of 100 
loads per load set). 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1 

O 

$NP'LDS 

or 

$PLDS 

AlO 

Keyword requesting the NPLDS or PLDS path of the LOAD 
module and introduces the NPLDS/PLDS card input. 

11-20 

print 
. or 

no-print 

AlO 

Visual keyword Indicating that a print key may appear in 
Col. 21-30, 31-40 or 41-50. If printing is desired a 
keyword must appear in the appropriate column (s). 

21-30 

MAP 


Keyword requesting the general print summary. 

31-40 

niiii 


Keyword requesting the general print summary and the 
output matrices (as placed on NPTAP or PTAP) . 

41-50 

CHEckOUt 

AlO 

Keyword requesting the general print summary, output 
matrices (as placed on NPTAP or PTAP) , and intermediate 
calculations. 
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Repeat cards 5. 1 through 5.6 for each required surface. 


Card 5.1 - Surface Definition 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

IrlO 

Surface 

AlO 

Keyword to identify the surface. 

11-15 

IS 

15 

Surface nujnber requested. 

(A default maximum number of 20. See Card 3.5 
for increase) 

16-20 


5X 

Not Used 

21-30 

CARDS 

or 

TAPE 

or 

blank 

AlO 

Keyword indicating the source of the mass data for 
surface IS. 

Keyword CARDS will cause the mass matrix to be read 
in Card 5.5. . 

Keyword TAPE will cause the mass matrix to be read from 
the tape name defined in the General Loads Input for 
MASSTP . 

The field is left blank if $PLDS is used on Card Set 5.0.- 


Card 5.2 — Comment (Optional) 

Repeat card 5.2 as desired to label the card data and printed output. 


mm 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-2 

C 

A2 

Keyword indicating a comment card that will be printed. 

3-70 


A8,6A10 

Available for comments 
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Card 5.3 — Surface Option Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

NODE 

or 

NODE -ALL 

AlO 

Keyword indicating which structural nodes of surface IS 
will be used to calculate PLDS or NPLDS. 

Keyword NODE means that specific nodes will be read 
on Card 5.4. 

Keyword NODE-ALL means all nodes defined in INTERP (L215) 
(Ref. 5) will be used. 

11-20 

OPTl* 

or 

OPT2** 

or 

blank 

or 

PLDS 

AlO 

Keyword defining the aerodynamic force interpolation method.. 
Keyword OPTl requires the P matrix to be read using 
Card 5.6. Program NPLDS is used. 

Keyword OPT2 requires the area of the structural panels to 
be read using Card 5.7. Program NPLDS is used. 

If the keyword is left blank, the structural and aero node 
locations are assumed to be the same. 

Keyword PLDS indicates that the structural and aero node 
locations are assumed to be the same and final results 
are written on the PLDS tape (without (M^ lias well as the 
NPLDS tape. 

If Card Set 5.0 is $PLDS, this keyword is left blank. 

21-30 

SCALRl 

ElO.O 

Scale factor to be applied to[R 2 l 
Default: SCALRl =1.0 

31-40 

SCALR2 

ElO. 0 

Scale factor to be applied to [M^l and iM^l 
Default : SCALR2 = 1.0 

41-50 

SCALR3 

ElO.O 

Scale factor to be applied to[$). 
Default; SCALR3 =1.0 

51-60 

PRINT 

or 

NO-print 

AlO 

A keyword to define the sequential print option after each 
surface calculation is performed. 

Keyword NO indicates no printing unless selected on 
$NPLDS/$PLDS (Card Set 5.0). 

Default 5 PRINT 


*OPTl inteipolates using a weighting factor matrix P. 

**OPT2 interpolates using a form of surface spline and requires the areas of the structural 
nodes to be input. 
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Card 5.4 — Structural Node Numbers (Optional) 

Omit card 5.4 if keyword on card 5.3 is NODE-ALL. 

Repeat card 5.4 as necessary until all nodes required are defined. 


COLS. 

KEyWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

6-10 

I INODE (1) 
I INODE (2) 

1415 

Numbers of structural nodes (monotonic-increasing) specifically 
requested to calculate PLDS and/or NPLDS loads. 


• 


For example; 1357 would include 1,3, 5, 7. If a string of 

• 

■ 

66-70 

• 

IINODE (14) 


nodes is required in addition to some selected nodes, it might 
be requested by 1 3 5 -10 13 14 where 5 -10 is interpreted 

to mean nodes 5, 6, 7, 8, 9 and 10. 


Card 5.5 — Diagonal Elements of Mass Matrix (Optional) 

This card is included only if keyword on card 5. 1 is CARDS. 
Repeat card 5.5 as necessary until all mass data required is used. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

m 

XMASS (1) 

7E10.0 

Diagonal mass matrix elements for surface IS corresponding 

. 

• 


to node numbers requested on Card 5.4 and/or Card 5.3. 

61-70 

a 

XMASS ( 7) 


Repeat until all elements are defined. 
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Card 5.6 — Weighting Factor Matrix P for Interpolation of Aero Forces to Structural Nodes* 
(Optional) 

This card is included if keyword on card 5.3 is OPTl *. 

The P-matrix (matrix relating aerodynamic nodes to structural load locations) is input as a 
sparse matrix. This card is repeated as necessary to include all rows and all columns per row 
required for the P-matrix (the P-matrix is an N x M matrix where N =the number of structural 
nodes and M = the number of aero panels on surface IS). 

Repeat card 5.6 as necessary until all P-matrix nonzero elements are read. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

IROW 

15 

Row number 

6-10 


5X 

Not Used 

11-15 

ICl 

15 

Column number 

16-20 


5X 

Not Used 


^1 


P-matrix element corresponding to column ICl**. 


IC2 

15 

Column number 

36-40 


5X 

Not Used 

41-50 

P2 

ElO.O 

P-matrix element corresponding to column IC2. 

51-55 

IC3 

15 

Column number 

56-60 


5X 

Not Used 

61-70 

P3 

ElO.O 

P-matrix element corresponding to column IC3. 


*See section 4.3.3, program NPLDS. 

**Py is a fraction of aerodynamic force at node j acting on structural node i. 
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Card 5.7 — Structural Node Areas for Interpolation of Aerodynamic Forces to Structural 
Nodes* (Optional) 

This card is included only if keyword on card 5.3 is OPT2*. 

Repeat card 5.7 as necessary until all areas of structural nodes requested are input. 


COLS . 

KEYWORDS/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

11-20 

21-30 

AREA (1) 
AREA (2) 
AREA (3) 

7E10.0 

Structural node areas corresponding to the node numbers 
requested, for surface IS. 

Repeat as necessary to include all areas for structural 
nodes defined by Cards 5.4 and/or 5.3. 

61-70 

AREA (7) 


♦Warning* This option is only valid over regions where 

i 2 

and are equal to or approximately equal to 

zero 


♦See section 4.3.3 program NPLDS. 


Card 5.8 - End of NPLDS/PLDS Card Input 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SEND 

AlO 

Keyword indicates the end of Card Set 5.0. 
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6.3.4 PROGRAM VBMT 


Card set 6.0 contains instructions and data required to form shears and bending moment 
loads using data read from the interpolation (L2I5) tape (SATAP), (ref. 5), and the 
equations of motion (L217) tape (EOMLOD), (ref. 4). 

Each of the load coefficient matrices cannot exceed a size of 100 x 70 (100 loads, 70 degrees 
of freedom). If more loads are required, the user must use more load set cards (card 6.3). 
The final output file has successive load sets separated by end-of-file marks. 

The order of input for each set 6.0 is displayed in figure 10. 

Omit card set 6.0 if no shear nor bending moment load coefficient matrices are required. 



i NO 
( END 1 


Figure 10.— Flow of Card Set 6.0 Input Data 
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Card. Set 6.0 — VBMT Load Option 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SVBHT 

AlO 

Keyword requesting the VBMT path of the LOADS module and 
introducing the VBMT card input. 

11-15 

print 

or 

no-print 

AIO 

Visual keyword indicating that a print key may appear in 
Col. 21-30, 31-40 or 41-50. If printing is desired a 
keyword must appear in the appropriate column (s). 

21-30 

MAP 

AlO 

Keyword requesting the general print summary. 

31-40 

MATrices 

AIO 

Keyword requesting the general print summary and the 
output matrices (as placed on LTAP) . 

41-50 

CHEckout 

AlO 

Keyword requesting the general print summary, output 
matrices (as placed on LTAP) , and intermediate calculations. 


Card 6.1 — Surface Sequence 

This sequence order of surfaces is used as a guide and a check on the order of card 6.4. The 
sequence order is variable except that if a dummy load defined and calculated on surface “n” 
is used in calculating loads on surface “m”, the surface sequence order must have surface n 
before m. This is necessary so that the dummy load is defined and calculated before being 
used in other load calculations. The required calculation sequence for surfaces is also the 
sequence for the final output tape. 

Repeat this card as necessary. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-20 

Surface 

sequence 

A10,10X 

Keyword to identify the surface sequence card. 

21-25 

• 

ISEQ (1) 

1015 

Surface numbers indicating the surface numbers to be 
processed and the proper sequence to calculate the 
desired loads. 

• 

• 

66-70 

ISEQ (10) 


(A default maximum number of 20 surfaces. See Card 3.5 
for increase) 
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Card 6.2 — Dummy Node Increase (Optional) 

This card increases the number of permitted dummy nodes. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

^xdum 

AlO 

Kayword requesting an increase in the maximum number of 




permitted- dummy nodes. 

11-15 

MXDUM 

15 

The number of dummy nodes permitted in this VBMT run. 




Defaul't • MXDUM * 20 


Repeat cards 6.3 through 6.9.4 for successive load sets. 
Card 6.3 — Load Set 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

LOAD- set 

AlO 

Keyword LOAD-set introduces the beginning of data cards 




for a new load-set. 

11-15 

ILS 

15 

Load-set number. 
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Repeat cards 6.4 through 6.9.4 as necessary to calculate loads on all surfaces. 


Card 6.4 — Surface Definition 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SURface 

AlO 

Keyword to introduce the surface number. 

11-15 

IS 

15 

Surface number requested. 

(A default maximum number of 20. See Card 3.5 for increase) 

21-30 

SCALE 1 

£10.0 

Scale factor applied to all structural node-summed shears 
and moments for^M^I (not applied to dummy nodes) . 

Default :SCALE1 - 1.0 

31-40 

SCALE2 

£10.0 

Scale factor applied to all aerodynamic node response force 
summed shears and moments for[M^>[Mg] (not applied to dummy 
nodes) . 

Default :SCALE2 =1.0 

41-50 

SCALES 

ElO.O 

Scale factor applied to all aerodynamic node gust force 
summed shears and moments for['|'] (not applied to dummy 
nodes) . 

Default :SCALE3 =1.0 

51-60 

PRINT 

or 

NO-print 

AlO 

Keyword to define the sequential print option after each | 

surface calculation is performed. 

Keyword NO indicates no printing unless selected on $VBMT 
(Card Set 6.0) 

Default : PRINT 
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Card 6.5 — Comment (Optional) 

Repeat card 6.5 as desired to label the card data and printed output. 


COLS. . 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-2 

C_ 

A2 

Keyword Indicating a comment card which will simply be 




printed. 

3-70 


A8,6A10 

Available for comments. 


Card 6.6 — Aero Panel Force Node Z Coordinate Overwrite (Optional). 

Aerodynamic slehder bodies and interference bodies are modeled as straight tubes. This 
option will replace the Z coordinates of the aero panel nodes from the equations of motion 
file, EOMLOD, with Z coordinates input on cards (normally, the actual structural locations) 
This does not change the magnitude or direction of the slender body and interference body 
forces. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

Z Coord 

■ 

AlO 

Keyword introducing the card to input the Z coordinates 
which will override selected aero-nodes received from EOM 
for surface IS. 

11-20 

ALL 

AlO 

Keyword indicating the z overwrite option for the nodes on 
surface IS. 

Keyword ALL indicates all nodes to be overwritten 
using Card 6.6.1 

Default (blank) indicates only some nodes to be overwritten 
using Card 6.6.2 
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Card 6.6.1 — Z Coordinates for All Aerodynamic Nodes 
This card is included only if 2nd keyword on card 6.6 is ALL. 
Repeat this card as required to define all nodes on surface IS. 


COLS. 

KEYWORD/ 

VARIABLE 

■ 1 

FORMAT 


DESCRIPTION 

1-10 

Z (1). 

7E10.0 

The overwrite Z 

coordinates (reference axis) for all 

61-70 

• 

• 

• 

Z (7) 


aerodynamic node 

numbers for surface IS. 


Card 6.6.2 — Z Coordinates for Specified Aerodynamic Nodes 

This card is included only if 2nd keyword on card 6.6 is blank. Repeat this card as 
required. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

11-20 

NODE (1)* 
Z (1) 

15 

ElO.O 

Aerodynamic node numbers 

The overwrite Z coordinate (reference axis) for 
aerodynamic node (1) . 

21-25 

• 

15 


31-40 

• 

ElO.O 


41-45 

NODE (3) 

15 


51-60 

Z (3) 

ElO.O 



*Note; Node (1) is the first aerodynamic node number to be overwritten. 



























Repeat cards 6.7 through 6.9 for all load numbers and dummy load numbers for surface IS. 
One or more of cards 6.9.1, 6.9.2, 6.9.3, and 6.9.4 is required with each card 6.7. 

Card 6.7 — Load Station/Dummy Station 

The load/dummy load numbers must increase monatonically beginning with 1 for the total 
VBMT run. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

or 

Dload 

A5 

Keyword Load indicates a load station. 

Keyword Dload indicates a dummy load station (no 
transformation, no output) (maximum number of 20 
dummy stations, see Card 6.2 for increase). 

6-10 

LNODE 

15 

Load / dummy load number. 

11-20 

AX * 

6E10.0 

BS 

load station/dummy station 
coordinates in reference axis system 

21-30 

AY 

BBL 

31-40 

AZ 

WL 

41-50 

ATHETX 

Thetax (rotation of load relative to reference axis system) 
(deg) - not required for dummy node 

51-60 

ATHETY 

Thetay (rotation of load relative to reference axis system) 
(deg) - not required for dummy node 

61-70 

ATHETZ 

Thetaz (rotation of load relative to reference axis system) 
(deg) - not required for dummy node 

(Default value for Thetax, Thetay, Thetaz 0 degrees) 


*Note: Columns 1 1-70 give the location and orientation of the load station. 

Angles of the load axis system 6^^, dy, and are measured from the reference 
axis to the load axis using the Eular rotation order defined in card 3.4. 
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Card 6.8 — Load Components 

Omit this card if keyword of card 6.7 is Dload. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

IKEY (1)* 

6(A2,8X) 

Keyword indicating the load components required after 




the transformation, for the final output. The keyword 

. 

• 


can be any IKEY (I) or combination 

• 

• 


of IKEY (I) s in order of VX, VY, VZ , MXX, MYY, MZZ. 

51-60 

IKEY (6) 




IKEY (Keyword) 

Description 

VX 

Shear, x direction (load axis) . 

VY 

Shear, y direction (load axis) . 

VZ 

Shear, z direction (load axis) . 

MXX 

Bending moment, about the x axis (load axis) . 

MYY 

Bending moment, a)x>ut the y axis (load axis) . 

MZZ 

Bending moment, about the z axis (load axis) . 


Card 6.9 — Load Panel Summation Direction 

Use cards 6.9.1 or 6.9.2, 6.9.3 and 6.9.4 or combinations of same. 

See figure 1 1 for load panel summation example. 

Card 6.9.1 — Load Panel Summation— All (Optional) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

ALL 

AlO 

Keyword indicating all nodes/surface IS are to be summed 
for LOAD or DLOAD number (LNODE) defined on Card 6.7. 
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Load at X|_, V(_, z^_ on surface 2 uses the relationships 
Y greater than yg (Card 6.9.4) 

+ 0.30 X aerodynamic forces at aeordynamic 

nodes (5-8) (Card 6.9.2) 

+ Dummy load at dummy node (1) 

to calculate the required loads. (Card 6.9.3) 


Figure 1 1.— VBMT Panel Summation Example 
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Card 6.9.2 — Load Panel Summation (Optional) (see fig. 11) 


COLS. 


1-10 


11-20 


21-30 

31-40 


41-50 


KEYWORD/ 

VARIABLE 


A.F. 


XYZ 


A.F. 

XYZ 


A.F. 

XYZ 


FORMAT 


A5.5X 


ElO.O 


A5, 5X 
ElO.O 


A5, 5X 
ElO.O 


DESCRIPTION 


Keyword indicating which nodes on surface IS will be 
included in the LOAD (or OLOAD) calculation for load 
LNODE (defined on Card 6.7). The periods are part of 
the keyword and the table below shows the possible 
choices for A and F. The A in the keyword indicates a 
direction along one of the axis in the reference axis 
system. The F in the keyword indicates which nodes in 
the A direction will be included in the load calculation*] 
(see example below) . 


The value of the x, y, or z coordinate (reference axis 
system) indicated by the A portion of the keyword above* 
(see example below) . 


A may be 

F may be 

X = X direction 

LT = less than 

y = y direction 

LE = less than or equal 

z = z direction 

EQ = equal 


GE = greater than or equal 


GT = qreater than 


Example: 

X.GT. 1 0.0 This means all nodes which have an x coordinate (reference axis) greater 
than 10.0 will be included in the load calculation. 

*Note; The user can limit the choice of nodes in the x, y, and z direction by using com- 
binations of X, y, and z for A. However, banding of coordinates is illegal, i.e.; 
Y.GT. and Y.LT. on same card 6.9.2. 
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Card 6.9.3 — Load Panel Summation— Specific Panels (Optional) 


Repeat this card as necessary. 



KEYWORD/ 



COLS. 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

STRucturep 

or 

AERop 

AlO 

Keyword defining Structural Panel or Aero Panel. 

11-20 

SCALE 

ElO.O 

A scale factor applied to the loads from nodes defined 
in Col. 21-60. 

Default 5 SCALE = 1.0 

21-25 

• 

• 

66-70 

IINODE (1) 

• 

• 

IINODE (10) 

1015 

Nodes/surface IS to be used in summing for shears and 
moments for LOAD or DLOAD number (LNODE) defined on 
Card 6.7. Each node required is input (monatonically 
increasing) . A negative number may be used to indicate 
a string of numbers; e.g., 1 -5 would mean node numbers 

1. 2, 3, 4, 5. 


Card 6.9.4 — Load Panel Summation— Dummy Nodes (Optional) 
Repeat this card if necessary. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

ADD dload 

AlO 

Keyword identifying a dummy load which must be included 
in the summation for LOAD or DLOAD number (LNODE) defined 
on Card 6.7. 

11-15 

LNODED (1) 

15, 

Dummy node number which must be included. 
Scale factor to be applied to the dummy load. 
Default : SCALED =1.0 

21-30 

SCALED (1) 

ElO.O 

31-35 

41-50 

LNODED (2) 
SCALED (2) 

15, 

ElO.O 

51-55 

61-70 

LNODED (3) 
SCALED (3) 

15, 

ElO.O 
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Note: If SCALED = +2.0 (symmetric case), the program sets Fyp, Ay = 0 and 

multiplies Exj), Fzd» ^YD where F = force, M = moment, D = dummy 

load, X, Y, Z is direction or axis. 

If SCALED = -2.0 (anti-symmetric case), the program sets Fxq, Fzq, Myp = 0 and 
multiplies FyQ, Mxp, and Mzp by 2.0. 

This will properly sum the left and right thin bodies for total fuselage loads for the 
symmetric or anti-symmetric analyses, i.e., sums left and right horizontal stabilizer 
onto aft body. 

Card 6.10 - End of VBMT Card Input 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SEND 

AlO 

Keyword indicates the end of Card Set 6.0. 


Card Set 7.0 — Job Terminator 
Always the last card of L218 (LOADS). 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SOUIT 

AlO 

Keyword indicating the last data for the LOADS module 
has been read (signals end of execution) . 


64 

















Bequirements or Function 


Input/Output File Names IFILE 


Introduce Hass Data 


Surface Card 


Hass Data 


Introduce Inertia Data 


Surface Card 


Key Words and/or Variables 


MASS 


Surface IS 


XHASS(l) XMASS(2) 


JMAT 


SURface IS 


XMASS(i) 


Reference 

Card Format Card Set CS 


SLOADS 

AlO I 

■ 

STITLE Title 

A10,6A10 

I : 

SGEN 

AlO 



6A10 


AlO 


A10,I5 


Inertia Data Type 
and Jormat 


Inertia Data for JTYPE 
FULl on CS 3.3.2 


Inertia Data For JTYPE 
SPArse on CS 3.3.2 


Transformation Order 


Increase Maximum 
Number of Surfaces 


XJMAT(l) XJMAT(2) 


NODE(l) XJMAT(l) 


TRAN sform lORDER 


HAXSUR ISMAX 


Non-Standard Buffer Size I LBUF LBUF 


End of General 
Card Input 


XJMAT(i) 


NODE(i) XJMAT(i) 


2 AlO 


7E10.0 


3(I5,5X,E10.0) 


2 AlO 


AlO, 15 


AlO, IS 


AlO 












































Requirements or Function 



Key Words 

and/or Variables 



Card Format 

Reference 
Card Set(CS) 

AVD Option 
Summary Print Option 

$AVD 

print 

no-print 

MAP 

MATrices 

CHEckout 



5A10 

4.S 

Surface Card 
Input Option 

Surface 

IS 

cards! 

TAPE . j 

I ROWS 

I COLS 



A10,I5,SX 

A10,2I5 

4.1 


C_ 

Comment 






A2,AB,6A10 

4.2 

Matrix Selection 

A-F 

(A=TA,TV 

,TD,RA,RV, 

RD and F'XYZ,XY,XZ,YZ 

X,Y,Z) 


6(A2,1X,A3,4X) 

4.3 

Surface Options 

NODE 

NODE-ALL 

LOCAL 1 
INERT 1 

SCALRl 

SCALR2 

SCALR3 

PRINT 

NO-print 


2A10,3E10.0, 

4.4 

COORD 

ANGLES 1 





AlO 


Structural Node Numbers 
NODE or. CS 4.4 

IINODE(l) IINODE(2) . . 

• 


IINODE(i) 

1415 

4.5 

Transformation Angles 
ANGLES on CS 4.4 

THETAX 

THETAV 

THETAZ 





3E10.0 

4.6 

Interpolation Option 
COORD on CS 4.4 

X Y 

z 

ALB 

ALT 

ALTT 

ACODE ISUB 

6E10.0,A5,I5 

4.7 

Modal Data 
CARDS on CS 4.1 

♦ (1) 

♦ (2) 

• • 

• 


♦ (i) 


7E10.0 

4.B 

End of AVD Card Input 

SEND 

AlO 

4.9 











































Requirements or File 


VBMT Option 
Summary Print Option 


Surface’ Sequence 


Increase Maximum 
Number of 
Dummy Nodes 


Surface Definition 


Key Words and/or Variables 


Z Overwrite 


Z Coordinates 
ALL on CS 6.6 


Z Coordinates 
blank on CS 6.6 


Load Station/ 
Dummy Station 


Load Components 


$VBMT 

print 

MAP 

MATrices 

CHEckout 








Surface 

sequence 

ISEQ(l) 

ISEQ(2) 

• • • 

ISEQ(i) 

I4axduin 

MxnuM 





LOAD- set 

ILS 





SURface 

IS 

SCALEl 

SCALE 2 

SCALES 

PRINT 

NO-print 


C_ Comment 


5Knk 


Z(l) Z(2) . . . Z(i) 


NODE(l) Z(l) N0DE(2) Z(2) ... NODE(i) Z(i) 


LNODE AX AY AZ ATHETX ATHETY ATHETZ 


IKEY(l) (VX,VY,VZ,HXX,HYY,MZZ) . . . IKEY(i) 


Card Format 


5A10 


A10,10X,10I5 


AlO.IS 


A10,IS 


A10,I5,SX, 

3E10.0,A10 


A2,A8,6A10 


2A10 


7B10.0 


3(I5,5X,E10.0) 


AS,IS,6E10.0 


6(A2,BX) 

















































6.4 MAGNETIC FILES (TAPE OR DISK) INPUT DATA 


6.4.1 EOMLOD 

L218 (LOADS) reads aerodynamic force matrices from the file EOMLOD which is generated 
by the program L217 (EOM) (ref. 4). L217 (EOM) writes the matrices on EOMLOD in the 
WRTETP format (see ref. 1) with one logical file per surface. Surface one is in file one, 
surface two in file two, etc. The structure of a typical file is displayed in figure 12. Please 
note that [Fp^ «g] is stored in the standard FORTRAN COMPLEX fashion. 


Size 


Matrix 


Description 


NNODES X 2 

NNODES X 1 
NNODES X NMODES 

NNODES X NMODES 

2 (NNOES) X NMODES 



( X and y local corrdinates of the 
surfaces nodes 

{Areas of surface elements 

{ Real part of 
aerodynamic forces 

( Imaginary part of 
aerodynamic forces 

( Gust forces 
(complex) 




I Repeat per 
\ reduced 
j frequency 


where; 


NNODES = Number of nodes on the surface < 100 
NMODES = Number of modes ^70 

Figure 12— Contents of a File on EOMLOD 


6.4.2 SATAP, 

L218 (LOADS) reads modal data from the file SATAP which is generated by the program 
L21 5 (INTERP) (ref. 5). L21 5 (INTERP) writes the modal data matrices onto SATAP in 
the WRTETP format (see ref. 1) with one logical file per surface. Surface one is in file 
one, surface two in file two, etc. The structure of a typical file is displayed in figure 13. 
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Size 

NNODES X NMODES 
NNODES X NMODES 
NNODES X NMODES 
NNODES X NMODES 
NNODES X NMODES 
NNODES X NMODES 
NNODEX X 6 

Variables x 1 


Matrix 



Description 

Modal translation 

Modal translation y 

Modal translation z 

Modal rotation about x-axis 

Modal rotation about y-axis 

Modal rotation about z-axis 

BS, BBL, WL, 6 ^, 6 , and 6 ^ for each node 
on the surface ^ 

Modal interpolation definition 
(see ref. 5 for description 
and length of the array) 


where; 

NNODES = Number of nodes on the surface < 1 00 

NMODES = Number of modes ^70 


Figure 13— Contents of a File on SATAP 


6.4.3 MASSTP 

The mass data for surface IS is found in logical file IS of MASSTP. Figure 14 displays the 
required matrix for a surface with N nodes. The program can select a subset of the N nodes 
via card input data. The matrix is in the READTP/WRTETP format described in reference 1. 
The mass data appears as a column matrix because only the diagonal elements are defined. 

Note: MASSTP may be the same as JMAT (sec. 6.4.4). The first matrix in each logical 
file of both is a mass matrix. 
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Size 


Matrix 


Description 


N X 1 



Masses of the N nodes on surface IS 


Figure 14.— Contents of a File on MASSTP 


6.4.4 JTAPE 

The inertia (J matrix) data for surface IS is found in logical file IS of JMAT. All matrices 
are in the READTP/WRTETP format described in reference 1. Figure 15 displays the 
matrices required for a surface with N nodes. All arrays must be present. The program can 
select a subset of the N nodes via card input. 

6.5 OUTPUT DATA 


6.5.1 PRINTED OUTPUT 

Normal sequential output is automatic unless suppressed in the submodule input (see card 
4.4 for AVD, card set 5.3 for NPLDS/PLDS, and card set 6.4 for VBMT). 

General output is a print summary selected by card set 4.0, 5.0, and 6.0 of each submodule, 
respectively, (AVD, NPLDS/PLDS, VBMT). MAP gives the general print summary for a 
given submodule. MATRICES will print, in addition to the general print summary, the 
matrices as they are written on the final output tapes. 

During job execution, the general output is written on the file lOUT. The last item L218 
does is copy lOUT to the regular output. In case of a program or system error, the contents 
of OUT should be copied to output (see sec. 6. 1 on control cards and notice the cards 
following EXIT). 

The printed output is illustrated in section 7.0 (sample problem). 
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Size 
NN X 1 

NN X 1 
NN X 1 
NN X 1 

NN X 1 

NN X 1 

NN X 1 

NN X 1 

NN X 1 


Matrix 



Description 

Mass data at structural nodes 

Mass static moments at structural 
nodes 


Rotary mass moments at inertia at 
structural nodes 


Mass products of inertia at 
structural nodes 


1 ^NN <N for a given matrix 

Figure 1 5.— Contents of a File on JTAPE 


6.5.2 MAGNETIC FILES (TAPE OR DISK) OUTPUT DATA 

L218 (LOADS) can generate different magnetic files. Each is described in one of the follow- 
ing subsections. All are written in the READTP/WRTETP format described in reference 1 . 

6.5.2. 1 AVDTAP 

The AVD loads path generates AVDTAP. There is one logical file for each load set. A 
typical logical file is displayed in figure 1 6. 

The physical size of the matrices Mj, M 2 and M 3 are equal and are appropriately loaded 
with zero rows as illustrated on the following page. 
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Size 


Matrix 


Description 


30 X 1 


NLD X NMODES 
NLD X NMODES 
NLD X NMODES 



DYLOFLEX header matrix described 
in table 1 

Displacement load coefficient matrix 
Rate load coefficient matrix 
Acceleration load coefficient matrix 


where: 



NLD 

= 

Number of loads ^ 100 

NMODES 

= 

Number of modes ^ 70 


Figure 16.— Contents of a File on A VDTAP 


In general there may be several surfaces in a load set. Then: 


[Mj] = 


0 


0 

0 


^surface 2 

Surface 1 


0 

0 


0 

0 


Mo 


surface 2 

Surface 2 

; [M2] = 

-0 

0 


0 



• 
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Table 1. — DYLOFLEX Header Matrix Contents 


Word 


Contents 


1 

2 

3 

4 

5 

6 
7 


7HYDL0FLX 

6HL218VC 

Run Date (YR/MO/DAY) 

Number of modes 
Number of load equations 
Number of gust penetration zones 

Number of frequencies at which aerodynamics are defined 


8 

0 

9 

0 

10 

0 

11 

0 

12 

0 

13 

0 

14 

0 

15 

0 

16 

0 

17 

0 

18 

0 

19 

0 

20 

0 

21 

i for 

22 

i for M 2 

23 

i for Mg 

24 

i for M^ 

25 

i for Mg 

26 

0 

27 

0 

28 

0 

29 

0 

30 

i for (j) 


i> 0 


i = 0 


indicates that the matrix 
is in the file 


indicates that the matrix 
is null and is not in the file 
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[M3] = 


M 

0 

0 

M: 

0 

0 

Me 


^surface 1 


surface 2 


surface 3 


6.5.2.2 NPTAP/PTAP 

The NPLDS loads path generates NPTAP. The PLDS loads path generates PTAP. NPTAP 
and PTAP are identical, except PTAP does not contain [M 3 ] . There is one logical file for 
each load set. A load set is the result of processing all surfaces requested following 
$NPLDS/$PLDS card, but before the next $ card. A typical logical file is displayed in 
figure 17. 

6.5. 2.3 LTAP 

The VBMT loads path generates LTAP. There is one logical file for each load set. A load 
set is the result of processing all surfaces requested following the LOAD-SET card, but 
before the next LOAD-SET or $ card. A typical logical file is displayed in figure 18. 

6.6 RESTRICTIONS 

1. The maximum size of any final output matrix is 100 x 70 

2. The maximum size of any input matrix is 100 x 70 

3. Units in L218 must be consistent with units from SATAP/EOMLOD 

4. The maximum number of reduced frequencies (k) is 20 

5. The maximum number of gust panels is 35 



Size 


Matrix 


Description 


30 X 1 

NLD X NMODES 
NLD X NMODES 

NLD X NMODES 

2 (NLD) X NGCP 



DYLOFLEX header matrix described 
in table 1 

Acceleration load coefficient matrix 

Load coefficient for ' 

generalized coordinate 

rate 

Load coefficient for 

generalized coordinate ' 

acceleration 

Load coefficient for 

excitation function 

force (complex) , 


where: 


NLD 

NMODES 
NGCP 
and (M 3 I 


Number of loads <100 
Number of modes <70 
Number of gust zones <35 
Surface 1 
Surface 2 
• 

Surface last 



similarly for [M^] , [Mg] and 


'Only on NPTAP 


Repeat for 
each 

frequency 
K=1, NK 


Figure 17.— Contents of a File on NPTAP/PTAP 



Size 


Matrix 


Description 


30 X 1 

NLP X NMODES 
NLDx NMODES 

NLD 

2(NLD)xNGP 

where: 

NLD 

NMODES 



DYLOFLEX header matrix described 
in table 1 


Acceleration load coefficient matrix 
(see fig. 19 for contents) 

Load coefficient for \ 

generalized coordinate I 

rate I 


Load coefficient for 
generalized 

coordinate acceleration 

Load coefficient for 
excitation function 
force (complex) 


f Repeat 
for each 
frequency 
K=1, NK 


Number of loads < 100 
Number of modes <70 
Number of gust zones <35 
Figure 18.— Contents of a File on LTAP 


NGCP 


[Wl o] 


’‘I 


pointy 




m 




point 


m 


> Surface 1 


M, 


point ^ 


M, 


><1 


M, 


m 


point 


m 


point ^ \ 


> Surface 2 


point ^ } 


Surface NS 


Similary for [M^l , [Mg] and [(p] at all frequencies. 


Figure 19.— Contents of on File LTAP 


6.7 DIAGNOSTICS 


The following diagnostics are printed by the program when errors are encountered. 
6.7.1 RGEN DIAGNOSTICS (GENERAL INPUT) 


Fatal Errors 

1 FETADD ERROR WHILE SETTING lOUT IN RGEN, IRR= 

2 THIS SHOULD HAVE BEEN A SLOADS CARD, BUT AMOD= 

3 FETADD ERROR - JTAPE, IRR= 

4 FETADD ERROR WHILE SETTING MASSTP IN RGEN, IRR= 

5 TRANSFORM CARD ILLEGAL, TXYZ= 

6 ERROR IN RGEN 1 , THE CARD READ WAS NOT A SURFACE CARD, AIS 

7 CALLED FROM RGEN 2 

THIS SHOULD HAVE BEEN A SURFACE CARD 

8 CALLED FROM RGEN 2 

J-MATRIX INPUT, SURFACE NO. LE. PREVIOUS NO.-IS, ISP= 

9 CALLED FROM RGEN 2 

TYPE OF J-MATRIX (SPARSE OR FULL) NOT SPECIFIED, AMOD, TYPE= 

10 CALLED FROM RGEN 2 

J-MATRIX, KEYWORD MASS, MX, IXX, ETC. BAD, AMOD= 

1 1 CALLED FROM RGEN 2 

J-MATRIX, IKEY CARD ERROR, IKEY, IKEYP= 

12 CALLED FROM RGEN 2 

J-MATRIX, FULL MATRIX SPECIFIED, NODES= 

1 3 CALLED FROM RGEN 2 
NODES EXCEEDS 100, NODES= 

14 CALLED FROM RGEN 2 
WRTETP ERROR, IRR= 



Warning Errors 


1 ERROR IN RGEN 1, ENCOUNTERED IN WRTETP,IR= 

2 CALLED FROM RGEN 2 

THE NO. ROWS IN THE J-MATRIX EXCEEDS NODES SPECIFIED NQDES, IROW= 
6.7.2 AVD DIAGNOSTICS 


AVD Program Diagnostics 


1 FATAL ERROR 

PROGRAM AVD ERROR NUMBER (IR) 
FETADD ERROR, IRR= 
where IR = 1 for file IFM 1 

2 for file IFM2 

3 for file IFM3 

4 for fUe IDISK 


AVDTAP Subroutine Diagnostics 


1 FATAL ERROR 

AVDTAP CAN NOT WRITE FINAL TAPE, FET IRR= 
(for file I AVD 1) 

2 WARNING MESSAGE 

AVDTAP DETECTS AN ERROR IN NUMBER OF 
ROWS FOR FINAL MATRIX 
12, ITOTAL= 


DISK Subroutine Diagnostics 

1 DYNAMIC STORAGE REQUIRES LOCATIONS 

DISK 1 Subroutine Diagnostics 
1 CALLED FROM DISK 1 
DISK 2 Subroutine Diagnostics 
1 CALLED FROM DISK 2 


DISKS 5 Subroutine Diagnostics 


Warning Message Errors 

I ERROR IN DISK WHILE READING MATRICES. 
LOCATION INDICATOR = IRR= 
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2 A NULL MATRIX WAS READ, THE LOCATION INDICATOR = 

3 THE NUMBER OF MODES READ IN MATRIX DOES NOT AGREE WITH PREVIOUS 

MODES, MXMODE = MXMOD = 

THE LOCATION INDICATOR = 

Note: IRR is the READTP error indicator 

MXMODE = current matrix 
MXMOD = previous matrix 
The location indicator defines; 


1 

— 

PHI-X matrix 

2 

- 

PHI-Y matrix 

3 

- 

PHI-Z matrix 

4 

- 

PHI-THETA-X matrix 

5 

- 

PHI-THETA-Y matrix 

6 

- 

PHI-THETA-Z matrix 

7 

- 

Geometry matrix 

8 

- 

SA matrix 

9 

- 

Geometry matrix 

10 

- 

Not used 

11 

- 

PHI-X matrix 

12 

- 

PHI-Y matrix 

13 

— 

PHI-Z matrix 

14 

- 

PHI-THETA-X matrix 

15 

— 

PHI-THETA-Y matrix 

16 

— 

PHI-THETA-Z matrix 

17 

- 

Geometry matrix 

18 

— 

Geometry matrix 
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MERGE Subroutine Diagnostics 


1 WARNING MESSAGE 

THE NUMBER OF ROWS IN THIS LOAD-SET EXCEEDS THE MAXIMUM 100 
PERMITTED* ITOTAL= 

MERGE 1 Subroutine Diagnostics 


(none) 


NTERP Subroutine Diagnostics 


1 WARNING MESSAGE 

ERROR CONDITION, INTER. GT. 6 
DETECTED IN SUB. NTERP 


NTERP 1 , NTERP 5 or NTERP 6 Subroutine Diagnostics 


(none) 

RAVD Subroutine Diagnostics 

1 DISK ERROR, CALLED FROM RAVD, IRR= 

RAVD 1 Subroutine Diagnostics 

1 TEMP 1 (I) = (keyword from input card) 

EAVDIA Subroutine Diagnostics 

1 ERROR IN RAVD I , TEMP = (keyword from input) 

RAVD 2 Subroutine Diagnostics 
(none) 

RAVD 3 Subroutine Diagnostics 

1 ERROR IN RAVD 3, NODES = 

NUMBER OF ANGLES = 

RAVD 4, RAVD 5, RAVD 5A, RAVD 5B, Diagnostics 
(none) 

RAVD 6 Subroutine Diagnostics 


(this is the error diagnostic subroutine for RAVD and associated routines) 


I 

1 ERROR NUMBER IN RAVD 

2 SURFACE CARD DID NOT HAVE S IN COLUMN 1 

3 ILLEGAL INFORMATION ON KEYWORD CARD 

4 END-OF-FILE ENCOUNTERED WHILE READING NODE CARDS 

5 END-OF-FILE ENCOUNTERED WHILE READING ANGLES CARDS 

6 END-OF-FILE ENCOUNTERED WHILE READING COORDINATE CARDS 

7 ILLEGAL OPTION (ANODE), NODE-AXIS CARD, CALLED FROM RAVD 

8 DISK READ ERROR, NODE-ALL OPTION, CALLED FROM RAVD 

9 RAVD lA TEST ON XYZ NOT SATISFIED 

10 RAVD 3, NODE INPUT INCONSISTENT WITH ANGULAR INPUT 
6.7.3 NPLDS/PLDS DIAGNOSTICS 

NPLDS Program Diagnostics 
1 FATAL ERROR 

SUBROUTINE NPLDS ERROR NUMBER 

FETADD ERRO, IRR = 

ERROR NUMBER 1 file lEOMLD 

2 file IDISK 

3 file NPTAP 

4 file MERGMB 
NPLDA Subroutine Diagnostics 

1 FATAL ERROR 

PROGRAM NPLDA, ERROR NUMBER 

FETADD ERROR, IRR= 

ERROR NUMBER 1 file MASSTP 

2 file IPTAP 

2 FATAL ERROR CALLED FROM NPLDA 

THE P-MATRIX WAS REQUESTED BUT NOT INPUT. 

GO TO THE NEXT INPUT. 
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3 READTP ERROR, CALLED FROM NPLDA, IRR =____ 

NPLDA 1 Subroutine Diagnostics 
(none) 

NPLDA 2 Subroutine Diagnostics 

1 CALLED FROM NPLDA 2, READING MASS TAPE, IRR = 

NPLDA 3, NPLDA 4 Subroutine Diagnostics 
(none) 

NPLDA 6 Subroutine Diagnostics 

(this is the error diagnostic subroutine for NPLDA and associated routine) 

1 ERROR NUMBER IN RAVD 

2 A SURFACE CARD DOES NOT HAVE AN S IN COLUMN 1 

3 ILLEGAL NODE OPTION (ANODE). 

4 ERROR READING DISK OR lEOMLD TO GET NODES. 

5 EOF ENCOUNTERED WHILE READING NODE CARDS. 

6 ERROR READING MASSTP TO GET MASS. 

7 EOR ENCOUNTERED WHILE READING P-MATRIX. 

8 IROW FROM MASSTP NOT EQUAL TO NODES. 

NPLDB Subroutine Diagnostics 

1 WARNING MESSAGE CALLED FROM NPLDB 

ERROR WHILE READING PZ IN NPLDB, IRR = 

2 WARNING MESSAGE CALLED FROM NPLDB 

ERROR WHILE READING GEOMETRY IN NPLDB, IRR = 

3 FATAL ERROR CALLED FROM NPLDB 

ERROR IN READING SATAP, IRR = 

4 FATAL ERROR CALLED FROM NPLDB 

ERROR IN READING EOMLOD, IRR = 



NPLDB 1 , NPLDB 2 Subroutines Diagnostics 


(none) 


NPLDD Subroutine Diagnostics 


WARNING MESSAGE CALLED FROM NPLDD 

1 ERROR IN READING LOCAL AERO COORDINATES FROM FILE 

lEOMLD, ICOL = 

2 ERROR IN READING LOCAL AREAS FROM lEOMLD, ICOL = ■ 

3 ERROR IN READING K-FREQUENCIES FROM lEOMLD, ICOL = 

4 THE NUMBER OF FREQUENCIES READ FROM lEOMLD EXCEEDS 

20, NK = 

5 SURFACE NUMBER = 

AERO AND STRUCTURAL NODES ARE ASSUMED IDENTICAL 

6 THIS OPTION (OPT 2) IS ONLY VALID OVER REGIONS WHERE DP/DX 
AND DP/DY ARE EQUAL OR APPROXIMATELY EQUAL TO CONSTATS. 

FATAL ERROR CALLED FROM NPLDD NODES = NAERO =. 

(when this diagnostic occurs it is in connection with number 5 above) 


NPLDD 1 Subrouting Diagnostics 


1 WARNING MESSAGE CALLED FROM NPLDP 1 ERROR IN DISK WHILE 

READING MATRICES. LOCATION INDICATOR = 

2 WARNING MESSAGE CALLED FROM NDLPP 1 A NULL MATRIX WAS READ, 

THE LOCATION INDICATOR = 


LOCATION INDICATOR 1 get local (x, y) coordinate 


2 get local area 


3 get frequencies 

4 ;;et (Fpl) matrix for M4 

5 get (fpl) matrix for M5 

6 get (Fpl) matrix for 0 

(location indicator is incremented for each (Fpl) 
matrix until matrices for all frequencies have been 
read) 



NPLDD 2 Subroutine Diagnostics 


(none) 

NPLDD 3 Subroutine Diagnostics 

There are 3 diagnostics, each preceded by 
WARNING MESSAGE CALLED FROM NPLDD 3 

1 TAPE READ ERROR ON IDISK WHILE READING SA FOR P-MATRIX 

CALC., IRR = 

2 PLATE 1 ERROR, INDS = 

3 NODES MUST BE .LE. NAERO. 

NODES, NAERO = 

NPLDD 4 Subroutine Diagnostics 

There are 2 diagnostics, each preceded by 

WARNING MESSAGE CALLED FROM NPLDD 4 

1 ERROR IN READING THE FORCE COEFFICIENT MATRIX. 

IROW, NAERO, K, I = 

2 ERROR IN READING THE FORCE COEFFICIENT MATRIX. 

ICOL, MXMODE, K, I = 

NPLDH Subroutine Diagnostics 

1 WARNING MESSAGE 

ERROR WRITING NPTAP, CALLED FROM NPLDH, IR, HL = 


NL 

Location 

File 

-1 

NPLDS 

PTAP 

-2 

PLDS 

PTAP 

1 

M 3 

NPTAP 

2 

M 4 

PTAP or NPTAP 

4 

1 

Is 

PTAP or NPTAP 

6 


PTAP or NPTAP 
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6.7.4 VBMT DIAGNOSTICS 


VBMT Program Diagnostics 
1 FATAL ERROR 

SUBROUTINE VBMT ERROR NUMBER 

SUBROUTINE VBMT ERROR NUMBER 

FETADD ERROR, IRR = 

ERROR NUMBER 0 file JTAPE 

1 file lEOMLD 

2 file IDISK 

3 file LTAP 

4 file MERGMB 

5 file INTMP 
VBMTA Subroutine Diagnostics 

Each of these diagnostics is preceded with 
FATAL ERROR CALLED FROM VBMTA 

1 AN EOF WAS ENCOUNTERED IN THE INPUT, CHECK YOUR INPUT 

2 AN S, C, Z, L OR D SHOULD HAVE BEEN IN COL I OF THIS CARD 

3 THE SURFACE NUMBER IS OUT OF SEQUENCE 

4 EXCEEDED MAX DUMMY NODES, ILD, MXDUM = 

VBMTA 1 Subroutine Diagnosti cs 

I FATAL ERROR CALLED FROM VBMTA I 

EOF ENCOUNTERED WHILE READING Z COORD 

VBMTA 2 Subroutine Diagnostic 

Each of these diagnostics is preceded with 
FATAL ERROR CALLED FROM VBMTA 2 

1 ERROR IN LOAD PANEL SUMMATION CARD, AIS = 

2 ONLY ONE (C. F. ) CARD ALLOWED 

87 


I 


3 


A (C. F.) CARD WAS READ, BUT NO NODES SELECTED 

4 (ADD DLOAD) WAS NOT DEFINED BY A DLOAD CARD 

5 LOAD SUMMATION INPUT ERROR, NODES, NAERO = 

VBMTA 3 Subroutine Diagnostics 

Each of these diagnostic is preceded with 
FATAL ERROR CALLED FROM VBMTA 3 

1 ERROR IN READING IDISK, IRR = 

2 ERROR IN READING SATAP, IRR ,= 

3 ERROR IN READING EOMLOD, IRR = 

VBMTC Subroutine Diagnostics 

(none) 

VBMTC 1 Subroutine Diagnostics 

Each of these diagnostics is preceded with 

FATAL ERROR CALLED FROM VBMTC 1 

1 ERROR IN READING JTAPE, IRR = 

2 ERROR IN READING JTAPE, ICOL, IROW = 

3 ERROR IN READING INTERP TAPE, IRR = 

4 ERROR IN READING INTERP TAPE 

IROW, ICOL, MXNODE, MXMODE = 

VBMTC 2 Subroutine Diagnostic 

1 FATAL ERROR CALLED FROM VBMTC 2 DUMMY NODE REQUESTED NOT 
DEFINED, NODED = 



VBMTC 3, VBMTC 4 Subroutine Diagnostics 


(none) 


VBMTD Subroutine Diagnostics 


1 WARNING MESSAGE CALLED FROM VBMTD ERROR IN READING K-FRE- 

QUENCIES FROM lEOMLD, ICOL = 

2 WARNING MESSAGE CALLED FROM VBMTD THE NUMBER OF FREQUENCIES 

READ FROM lEOMLD EXCEEDS 20, NK = 


VBMTD 1 Subroutine Diagnostics 


1 WARNING MESSAGE CALLED FROM VBMTD 1 ERROR IN DISK WHILE READ- 
ING MATRICES. LOCATION INDICATOR = IRR = 

2 WARNING MESSAGE CALLED FROM VBMTD 1 A NULL MATRIX WAS READ, 

THE LOCATION INDICATOR = 

LOCATION INDICATOR 1 read frequencies 

4 read(FpI) matrix for M4 

5 read (Fpl) matrix for Mg 

6 read (Fpl) matrix for 

(the location indicator is incremented for each (Fpl) 
matrix until matrices for all frequencies have been read) 

VBMTD 2 Subroutine Diagnostic 


(none) 


VBMTD 3 Subroutine Diagnostics 


1 WARNING MESSAGE CALLED FROM VBMTD 3 ERROR IN READING THE 

FORCE COEFFICIENT MATRIX. IROW, NAERO, I, I = 

2 WARNING MESSAGE CALLED FROM VBMTD 3 ERROR READING THE FORCE 

COEFFICIENT MATRIX. ICOL, MXMODE, K, I = 


VBMTD 4 Subroutine Diagnostics 


1 FATAL ERROR CALLED FROM VBMTD 4 DUMMY NODE REQUESTED NOT 
DEFINED, NODED = 
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VBMTF Subroutine Diagnostics 


1 WARNING MESSAGE 

ERROR WRITING LTAP, CALLED FROM VBMTF, IR, NL 
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7.0 SAMPLE PROBLEM 


In order to exercise most of the options in this program, a very simple, fictitious airplane 
was modeled. The airplane modeled was a wing mounted two engined T-tail transport 
airplane. The airplane is shown in figure 20 with the lumped masses and aerodynamic grid 
superimposed. The inertia data, structural node numbers, and coordinate locations are 
given in table 2; the mode shapes for the corresponding node numbers in table 3; and the 
aerodynamic node number, coordinate locations, ard forces in table 4. The structural data 
magnetic file was generated with L215 (INTERP). The aerodynamic data was written onto 
a magnetic file in the format normally generated by L2 1 7 (EOM). 

The inertia data was read into L218 (LOADS) by the RGEN module. AVD, PLDS, NPLDS, 
and VBMT load coefficient matrices were generated on all applicable surfaces by exercising 
their respective modules for a vertical gust analysis. The analysis consisted of two rigid body 
degrees of freedom (vertical translation and pitch) and one elastic degree of freedom. The 
printed output for all of the load matrices was edited to reduced the amount of pages, but 
yet to yield a representative sample. 

Boeing Commercial Airplane Company 
P.O. Box 3707 
Seattle, Washington 98124 
May 1977 
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Figure 20— Airplane 


Surface no. 



S3 




























'•O 

4^ 


Table 2. — (Concluded) 


Surface 3* (total fuselage) 



Structural Nodes 

8^ 

w 

Wx 

Wz 

Ixx 

lyy 

Izz 

Node No. 

B.S. 

-B . B . L . 

W.L. 

degs . 

lbs 

in lbs 

in lbs 

Ib-in^ 

Ib-in^ 

Ib-in^ 

1 

250 

0 

200 

mu 

50 

0 

-12000 

KBBHi 



2 

400 

0 

200 

■■ 

200 

0 

0 



K 

3 

600 

0 

200 

0. 

600 

0 

1200 

1X1 o' 

6X10' 

6X10' 

4 

850 

0 

200 

Km 

1000 

0 

-20000 

2X10' 

1X10* 

1X10* 

5 

1250 

0 

200 

KH 

1000 

0 

0 

2X10' 

1X10* 

1X10* 

6 

1500 

0 

200 


1000 

0 

0 

2X10' 

ixio* 

1X10* 

7 

1800 

0 

200 

0. 

700 

0 

-21000 

2X10' 

7X10' 

7X10' 

8 

2050 

0 

210 


600 

0 

6000 

1X10' 

6X10' 

6X10' 

9 

2250 

0 

230 


400 

0 

0 

1X10' 

4X10' 

4X10' 

10 (included 
vert, fin) 

2450 

0 

250 


300 

0 

15000 

2X10' 

2.5X10' 

2.5X10' 


Surface 4« 


1 

2450 

50 

400 

0. 

70 

mmm 

0 


1X10' 

Q SB 

2 

2460 

150 

400 


50 

HI 

■■ 

iT~^l 

1X10' 


3 

2470 

300 

400 


50 

1500 


1 ■ 

1X10' 


4 

2480 

450 

400 

BH 

20 

-600 

■■ 

2X10’ 

.4X10' 



Surface 5* 

Not required for vertical gust analysis. 


'See Figure 20 






















Table 3. - 


Surface 1* 



See Figure 20 


Mode Shapes 




>5 


















































Table 4. — Aerodynamics — Node Locations and Forces 


Surface 1* 





k = 

0 



1 

Aero 
Node NO. 

B.S. 

B.B.L. 

W.L. 

Fy(q) 

Fy(q) 

Fz (q) 

Fz(q) 

1 

950 

500 

ISO 

0 

0 5. 


0 

f !f 

mm 


2. .1 

2 


500 

180 

0 

0 2. 

1— IH 

0 



1. .05 








k = 

1 






1 

mmm 

500 


0 

0 

4. 

B 


.1 

B 

150 

10 

15 

1. 

.05 

2 


500 


0 

0 

2. 

D 

Hi 

.05 

0 

100 

4. 

5. 

.5 

• .01 







k = 

- 0 


k = 1 1 





Fy(g 

ust) 

Fz (qust) 

Fy(c 

lust) 

Fz (gust) 





R 

I 

R 

1 

R 

I 

R 

^ 1 

1 

mam 

500 


0 

Di 

■1 


0 

0 

15 

■Hi 

2 

BBI 

500 


0 


Hal 

■■ 

0 

0 

5 



See Figure 20 

'These aerodynamic forces are for the three degrees of freedoms 



















































(Continued) 















Table 4 (Continued) 


Surface 3 * 




Aero 


■H 








Fz(qust) 

Node NO. 

B.S. 


W.L. 


F(q) 



F(q) 


R 

I 

1 



200 

0 

3000 

-100 

3 

-20 

1 

3 


2 



200 

0 

200 

-20 

.2 

-2 

.1 

.2 


3 



200 

0 

10 

-10 

.01 

0 

mm 

.01 


4 

B^l 

0 

200 

0 

-30 

-10 

-.3 


n 

-.3 


5 

B^h 

0 

200 

0 

-20 

-5 

-.2 

.s 

0 

-.2 


6 


0 

200 

0 

70 

0 

.7 

1 

.1 

.7 

IDI 






0 

200 

.01 

300 

-100 

.3 

-2 

n 

.3 

.3 



0 

200 

.01 


-20 

.2 

-.2 


.2 

.2 



0 

200 

.01 

10 

-10 

.01 

0 


.01 

.01 




200 

.01 

-3 

-5 

-.003 

0 

D 

-.003 

-.003 

5 



200 

.01 

-2 

0 

-.002 

.05 

mm 

-.002 

-.002 


hSBI 


200 

.01 

7 

10 

.007 

.1 

H 

.007 

.007 


See Figure 20 





















Table 4. 



See Figure 20 























$LO*OS 

$T1TLE 

*T1TL£ VERIFICATION TEST CASES FOR L21f. 
ATITLE 


tCEN 

SATAP 

SATAPl 

EOHLUO 

EOHLOl 

HASSTP 

JNAT 

SURFACE 

1 

NODES 

1 


HASS 

FULL 




1000. 

IXX 

100000. 

FULL 




IVV 

200000. 

FULL 




Ill 

AOOOOO. 

FULL 




SURFACE 

2 

NODES 

12 


MASS 

FULL 




100. 

100. 

200. 

250. 

200. 

ASO. 

600. 

700. 

1000. 

1200. 

SURFACE 

3 

NODES 

10 


MASS 

FULL 




SO. 

200. 

600. 

1000. 

1000. 

600. 

600. 

300. 



M 

FULL 




-12000. 

0. 

12000. 

-20000. 

0. 

6000. 

0. 

15000. 



IXX 

FULL 




500. 

1000. 

10000. 

20000. 

20000. 

10000. 

lOOCO. 

20000. 



lYY 

FULL 




5000. 

20000. 

6CC00. 

100000. 

lOOCOO. 

60000. 

60000. 

25000. 



IZZ 

FULL 




5000. 

.20000. 

60000. 

100000. 

100000. 

60000. 

60000. 

25000. 



SURFACE 

6 

NODES 

6 


HASS 

FULL 




70. 

SO. 

SO. 

20. 


MX 

FULL 




-3500. 

-500. 

ISCO. 

-600. 


IXX 

FULL 




7000. 

5000. 

5000. 

2000. 


lYY 

FULL 




10000. 

10000. 

10000. 

6000. 


IZZ 

FULL 




10003. 

10000. 

lOCOO. 

6000. 


SURFACE 

5 

NODES 

3 


HASS 

SPARSE 




1 

100. 

2 

SO. 

3 

MX 

FULL 




-5000. 

1250. 

0. 



IXX 

FULL 




1000. 

500. 

500. 



lYY 

FULL 




2000. 

1000. 

1000. 






V6RT 

2 



VERT 

3 



VERT 

6 



VERT 

5 



VERT 

6 

JNAT 


VFRT 

7 



VERT 

8 



VERT 

9 



VERT 

10 



VERT 

11 



VERT 

12 



VERT 

13 



VERT 

16 



VERT 

15 



VERT 

16 



VERT 

17 



VERT 

18 



VERT 

19 

300. 

600. 

VERT 

20 



VERT 

21 



VERT 

22 



VERT 

23 

1000. 

TOO. 

VERT 

26 



VERT 

25 



VERT 

26 

0. 

-21000. 

VERT 

27 



VERT 

28 



VERT 

29 

20000. 

20000. 

VERT 

30 



VERT 

31 



VERT 

32 

100000. 

70000. 

VERT 

33 



VERT 

36 



VFRT 

35 

100000. 

70000. 

VERT 

36 



VERT 

37 



VERT 

38 



VERT 

39 



VERT 

60 



VERT 

61 



VERT 

62 



VERT 

63 



VERT 

66 



VERT 

65 



VERT 

66 



VFRT 

67 



VERT 

68 



VERT 

69 



VERT 

50 

SO. 


VERT 

51 



VERT 

52 



VERT 

53 



VFRT 

56 



VERT 

55 



VERT 

56 



VERT 

57 


uz 

FULL 







VERT 

58 

2000. 

1000. 

1000. 






VERT 

59 

tENO 

GEN 







VERT 

60 

$AV0 

HAP 

HATRICES 






VERT 

61 

SURFACE 

2 

CAROS 

2 

3 




VERT 

62 

TA-2 

TV-2 

TO-2 






VERT 

63 

COORD 


.00259 

1. 





VERT 

64 

1250. 

1000. 

200. 






VERT 

65 

1350. 

1000. 

20C. 






VERT 

66 

1. 

50. 

1. 






VERT 

67 

1. 

150. 

.8 






VERT 

68 

SURFACE 

1 

CARDS 

1 

3 




VERT 

69 

I A-X 

TV-X 







VERT 

70 

COORO 

LOCAL 

1. 

1. 





VERT 

71 

1000. 

500. 

18C. 






VERT 

72 

0. 

10. 

. 1 






VERT 

73 

SURFACE 

1 

CAROS 

1 

3 




VERT 

74 

TV-V 

lo-y 







VERT 

75 

COORO 



1. 





VERT- 

76 

1000. 

500. 

180. 






VERT 

77 

0. 

0. 

.1 






VERT 

78 

SURFACE 

1 

CAROS 

1 

3 




VERT 

79 

IA-2 








VERT 

80 

COORO 


.00259 






VERT 

81 

1000. 

500. 

180. 






VERT 

82 

1. 

-200. 

.4 






VERT 

83 

SURFACE 

1 

CAROS 

1 

3 




VERT 

84 

RV-X 








VERT 

85 

COORO 




1. 




VERT 

86 

1000. 

500. 

180. 






VERT 

87 

0. 

0. 

.01 






VERT 

88 

SURFACE 

1 

CAROS 

1 

3 




VERT 

89 

RA-Y 

RO-Y 







VERT 

90 

COORO 








VERT 

91 

1000. 

500. 

1(0. 






VERT 

92 

0. 

1. 

.02 






VERT 

93 

SURFACE 

1 

CAROS 

1 

3 




VERT 

94 

RA-2 

RV-2 

RO-Z 






VERT 

95 

COORD 




1. 




VERT 

96 

1000. 

500. 

180. 






VERT 

97 

0. 

0. 

.001 






VERT 

98 

SEND 

GENERAL 







VERT 

99 

$AVO 

HAP 

HATRICES 






VERT 

100 

surface 

1 

TAPE 






VERT 

101 

TA-XV2 

RV-XV2 







VERT 

102 

COORO 








VERT 

103 

900. 

500. 

280. 




PAR 


VERT 

104 

SURFACE 

2 

TAPE 






VERT 

105 

TA-1 

TD-2 







VERT 

106 

NODE-ALL 

LOCAL 







VERT 

107 

SURFACE 

3 

TAPE 






VERT 

108 

TA-2 








VERT 

109 

COORO 








VERT 

110 

ISO. 

0. 

250. 

-100. 

0.- 

50. 


1 

VERT 

111 

1000. 

0. 

200. 

ISO. 

0. 

0. 


9 

VERT 

112 

SURFACE 

3 

TAPE 






VERT 

113 


104 


rw 

RW-Y 






VERT 

114 

NODE 

LOCAL 






VERT 

115 

1 

3 4 

7 -10 





VERT 

116 

SURFACE 

4 

TAPE 





VERT 

IIT 

TA-Z 

RV-KV 

RO-X 





VERT 

118 

NODE 

ANGLES 






VERT 

119 

1 

3 






VERT 

120 

0. 

0. 

-4. 





VERT 

121 

0. 

0. 

-4. 





VERT 

122 

SEND 

AVO 






VERT 

123 

fAVO 

HAP 

HATRICES 





VERT 

124 

SURFACE 

2 

TAPE 





VERT 

125 

TA-Z 







VERT 

126 

COORD 

LOCAL 






VERT 

127 

1300. 

1000. 

20C. 





VERT 

128 

1350. 

900. 

20C. 





VERT 

129 

1000. 

100. 

200. 





VERT 

130 

900 ^ 

150. 

200. 





VERT 

131 

1000. 

190. 

200. 





VERT 

132 

SEND 

AVO 






VERT 

133 

tPLOS 

MAP 

HATRICES 





VERT 

134 

SURFACE 

2 

TAPE 





VERT 

139 

NODE-ALL 





PRINT. 


VERT 

136 

SURFACE 

4 

TAPE 





VERT 

137 

NODE 





PRINT 


VERT 

138 

1 

-3 5 

7 -9 





VERT 

139 

SURFACE 

1 

TAPE 





VERT 

140 

NOOE-ALL 





PRINT 


VERT 

141 

SEND 

AVO 






VERT 

142 

INPLOS 

HAP 

HATRICES 





VERT 

143 

SURFACE 

2 

TAPE 





VERT 

144 

NOOE-ALL 

OPT 2 




PRINT 


VERT 

149 

330. 

335. 

330. 

335. 

330. 

335. 

330. 

VERT 

146 

335. 

330. 

335. 

330. 

335. 



VERT 

147 

SURFACE 

4 

TAPE 





VERT 

148 

NODE 

OPTl 




PRINT 


VERT 

149 

1 

-3 






VERT 

150 

1 

4 

.5 

9 

.5 

6 

.5 

VERT 

151 

1 

7 

1. 

a 

1. 

9 

1. 

VERT 

152 

2 

1 

.5 

2 

.5 

3 

.5 

VERT 

153 

2 

4 

.5 

9 

.5 

6 

.9 

VERT 

154 

3 

1 

.9 

2 

.9 

3 

.9 

VERT 

155 

(END 

NPLOS 






VERT 

156 

ATITLE 

VBNT TEST CASE 





VERT 

157 

AVBNT 

HAP 

HATRICES 





VERT 

158 

SURFACE 

SEQUENCE 

1 

2 4 

3 



VERT 

159 

LOAD-SET 

1 






VERT 

160 

SURFACE 

1 

.00259 





VERT 

161 

C VERTICAL GUST 

747 AIRPLANE 

CHECK 

CASE NO.l 



VERT 

162 

C NACELLE 






VERT 

163 

LGAO 

1 1100. 

50C. 

200. 




VERT 

164 

VX 

VY 

VZ 

HXX 

HYY 

HZZ 


VERT 

165 

ALL 







VERT 

166 

OLOAO 

11200. 

900. 

200. 




VERT 

167 

ALL 







VERT 

168 

SUPFACE 

2 

.00259 





VERT 

169 



C HINC 







VERT 

170 

LOAD 

Z1200. 

900. 

200. 

5. 

0. 

-20. 

VERT 

171 

VX 

VV 

VZ 

NXX 

NVV 

NZZ 


VERT 

172 

V.GT. 

500; 






VERT 

173 

LOAD 

3 1050. 

190. 

200. 




VERT 

174 

VX 

VV 

VZ 

NXX 

NVV 

NZZ 


VERT 

175 

V.GT. 

ZOO. 






VERT 

176 

STAUCTUAEP .3 

11 

12 




VERT 

177 

AEKOP 

.9 

IT 

-20 




VERT 

178 

ADO CLOAO 1 






VERT 

179 

SURFACE 

4 

.0C259 





VERT 

180 

C HORIZONTAL STAGILIZER 





VERT 

181 

LOAD 

A 2430. 

0. 

400. 




VERT 

182 

VX 

VV 

VZ 

NXX 

NVV 

NZZ 


VERT 

183 

ALL 







VERT 

184 

OLOAO 

2 2430. 

0. 

400. 




VERT 

185 

ALL 







VERT 

186 

SURFACE 

3 

.00299 





VERT 

187 

C BODY 







VERT 

188 

Z COORO 







VERT 

189 

6 

240. 






VERT 

190 

LOAD 

5 800. 

0. 

200. 




VERT 

191 

VX 

VZ 

HVV 

NZZ 




VERT 

192 

X.LT. 

800. 






VERT 

193 

LOAD 

4 1800. 

0. 

200. 




VERT 

194 

VV 

VZ 

MXX 

NVV 




VERT 

195 

X.CT. 

1800. 






VERT 

196 

STRUCTUREP .9 

7 





VERT 

197 

AEROP 

.5 

9 





VERT 

198 

ADD CLOAO 

1 2 

2. 





VERT 

199 

LOAD 

7 1200. 

0. 

200. 




VERT 

200 

VZ 

HYV 






VERT 

201 

X.GT 

1200. 






VERT 

202 

AOO OLOAO 

2 

2. 





VERT 

203 

AENO 

V8MT 






VERT 

204 

tOUIT 







VERT 

209 


106 



• MOCItAM L2ISA1 VERStON APR 20t TT NOM RUANING. 

• Tl-E PROGRAM IS PART OF THE OVLOFLX SYSTEM 

• DEVELOPED FOR NASA UNDER CONTRACT NASl-lSGIR. 

• CATE OF RUN IS T1/OA/21. 

• TINE OF RUN IS IR.29.S2. 


( ILOADS 
I (TITLE 

I (TITLE VERIFICATION TEST CASES FOR L21(. 

< (TITLE 
I (GEN 

ISATAF SATAPl EONLOO EOMIDI HASSTP 


JHAT 


TAPE FUE NAMES hERE READ. THOSE NAMES TO BE USED ARE 

AVOTAF NFTAP PTAP LTAP EONLOl . SATAPt JNAT JTAPE 


IJMAT 


(SURFACE 1 

MOOES 

1 



J-NATRIX DATA 

FDR SURFACE 

1 ARE BEING READ 

(NASS 

FULL 





.1000EA04 




(UX 

FULL 





.lOOOEADA 




IIW 

FULL 





.2000EAOA 




1122 

FULL 





.AOOOEAO* 




(SURFACE 2 

NODES 

12 



J-MATRIX DATA 

FOR SURFACE 

2 ARE BEING READ 

(MASS 

FULL 





.1000EA03 

.lOCOE^OS 

.2000EP03 

.2900E»03 


.1000E04 

. ISOOE+OA 



(SURFACE 3 

NODES 

10 



J-NATRIX DATA 

FOR SURFACE 

3 ARE BEING READ 

(MASS 

FULL 





.SOOOEAD2 

.2000E*03 

.NOOOEi’OS 

.1000E»04 

(M2 

FULL 





-.I200E«05 

0. 

.1200E+0S 

~.2000E«0S 

IIXX 

FULL 





.S000EA03 

.1000E404 

.1000E»0$ 

.2000E»0S 

(IVY 

FULL 





. 5000E«04 

.20COE40S 

.AOOOE»09 

. 1000E*06 

(122 

FULL 





.9000E»04 

.2000E40S 

.M00E«0S 

. 1000E«OA 

(SURFACE A 

NODES 

4 



J-NATRIX DATA 

FOR SURFACE 

4 ARE BEING READ 

(NASS 

FULL 





.TOOOE«02 

.9OC0E402 

.M00E*02 

.2000E«02 


.2Q00E403 

.3000E»03 

.4000E»03 

1 

.4500E»03 

.6000E403 

.1000E*04 

.lOOOE+04 

t 

.7000E*03 

.6000E+03 

.400CE»03 

0. 

0. 

-.2I00E+05 

I 

.6000E*04 

0. 

.2000E«OS 

.2000E«OS 

1 

.2000E«0S 

1 

. 1000E«05 

.1000E40S 

.lOOOEfOA 

.lOOOE«OA 

.T000E*05 

1 

.6000EYOS 

.4000E405 

.lOOOEfOA 

.IOOOEaOA 

.TOOOE«OS 

.MOOE40S 

.4000E«0S 


.T000EA03 

.3000E»03 

.1900E«09 

.2000E»0S 

.2S00EP09 

.2500E*03 



(NX 

FULL 





-.3500E»04 

-.SOC0E*03 

,1S00E*C4 

-.6000E*03 

i IXX 

FULL 

.70006*04 

.SOCOE*04 

.M00E*04 

.2000E*04 

1 IVY 

FULL 

,1000E*0S 

.10C0E*CS 

.1000E*05 

.4000E*04 

1 IZZ 

FULL 

.1000E*0S 

.lOC0E*0S 

,IOOOE*CS 

.4000E*04 

(SURFACE S 

NODES 

3 



J-MATRIX DATA 

FOR SURFACE 

9 ARE BEING READ 

(MASS 

sparse 

.1000E*03 

•SOC0E*02 

.5000E*02 


(NX 

FULL 





-.5000E*04 

.1230E*04 

0. 


( IXX 

FULL 

. 1000E*04 

.90C0E*03 

•S000E*03 


( lYY 

FULL 

.2000E*04 

.10C0E*04 

.1000E*04 


I Ul 

FULL 

.2000E*04 

.10COE*C4 

.1000E*04 


1 SEND 

CEN 





THE J-HATRIX HAS BEEN INPUT BY CARDS 


LOM-SET 


1 


SURFACE 

NODE 

REF. 

COORDINATES ICHI 




RCH 

NO. 

NO. 

■S/X 

aSL/Y 

kL/2 

Ml 

M2 

N3 

NO. 

2 

1 

12SO.O 

ICOO.O 

200.0' 



T2 

1. 


2 

13SO.O 

ICOO.O 

2C0.0 



12 

2. 

2 

1 

12SO.O 

ICOO.O 

2C0.0 


T2 


3. 


2 

ISSO.O 

ICOO.O 

200.0 


T2 


4. 

2 

1 

1290.0 

ICOO.O 

2C0.0 

T2 



9. 


2 

1990.0 

ICOO.O 

200.0 

T2 



6. 


SURFACE 

NODE 

REF. 

COORDINATES ICNI 



RCH 

NO. 

NO. 

■S/X 

aai/v 

M./2 

HI R2 

H3 

NO. 

1 

1 

1000.0 

900.0 

100.0 


TX 

7. 

1 

1 

1000.0 

900.0 

100.0 

IX 


0 . 


SURFACE 

NODE 

REF. 

COORDINATES ICHI 



RCH 

NO. 

NO. 

■S/X 

OOL/V 

HL/2 

HI 

H2 

H3 NO. 

1 

1 

1000.0 

900.0 

lOO.O 


TV 


1 

1 

1000.0 

900.0 

100.0 

TV 


10. 


SURFACE 

NO. 

NODE 

NO. 

REF. 

■S/X 

COORDINATES ICNI 
■OL /V NL/2 

HI 

H2 

H3 

RCH 

NO. 

1 

1 

1000.0 

900.0 

100.0 



T2 

11. 

SURFACE 

NO. 

NODE 

NO. 

REF. 

os/x 

COORDINATES ICHI 
■OL/V HL/2 

HI 

H2 

N3 

ROH 

NO. 

1 

1 

1000.0 

900.0 

100.0 


RX 


12. 



SURFACE NODE 

REF. 

COORDINATES tCH| 




ROW 

NO. 

NO. 

BS/X 

BBL/V 

t|L/2 

HI 

N2 

H3 

NO. 

i 

i 

iOOO.O 

soo.o 

leo.o 



RV 

13. 

i 

i 

iOOO.C 

soo.o 

110.0 

RV 



14. 

SURFACE NODE 

REF. 

COORDINATES ICHI 




RCk 

NO. 

NO. 

BS/X 

BBL/V 

NL/2 

HI 

H2 

H3 

NO. 

i 

1 

1000. 0 

soo.o 

180.0 



RZ 

IS. 

i 

i 

1000. 0 

500.0 

180.0 


RZ 


16. 

1 

1 

1000.0 

soo.o 

180.0 

RZ 



IT. 

Hi 1 

0. 

0 


0. 





Hi 2 

0. 

0 

. 

C. 





HI 3 

0. 

0 

. 

C. 





Hi 4 

0. 

0 

. 

0. 






HI 

HI 

HI 

HI 

HI 

Hi 

HI 

Hi 

HI 

Hi 

Hi 

Hi 

Hi 


5 

6 
7 

e 0. 
*> 0 . 
10 0 . 
11 0 . 
12 0 . 

13 0. 

14 0. 

15 0. 

16 0. 
17 0. 


.lOOOOE+01 
. 10000E60 1 
0 . 


.50000E*02 

.i:OOOE»03 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

.1COOOE401 

0 . 

0 . 

0 . 


.10000E«01 

.80000E4’00 

0 . 

0 . 

0 . 

.lOOOCEtOO 

0 . 

0 . 

0 . 

.2O0O0E>0i 

0 . 

0 . 

.100006-02 


H2 

H2 

H2 

H2 

H2 

H2 

H2 

H2 


0 . 

0 . 

.lOOOOE^Ol 

.iOOOOE^Ol 

0 . 

0 . 

0 . 

0 . 


0 . 

0 . 

.5000CE+02 

.1SOOOE«03 

0 . 

0 . 

0 . 

.iC0O0E«02 


0 . 

C. 

.lOOOOE + 01 

.eooocE«^oo 

0 . 

0 . 

0 . 

.1COOCE«00 


lOOOOE^OO 


H2 9 0. <>• 

M2 10 0. 0. 

M2 11 0. 0. 

M2 12 0. 0. 

M2 13 0. 0* 

M2 U 0. 0. 

M2 19 0. 0. 

M2 16 0. 0. 

M2 IT 0. 0. 


ol 

c. 

ol 

c. 

c. 

o! 


lOOOOE-01 

lOOOCE-02 


MS 

1 

.29900E-02 

.129SCE400 

.2S9CCE-02 

M3 

2 

.29900E-02 

.3aasoE*oo 

.20720E>02 

M3 

3 

a. 

0. 

0. 

M3 

6 

0. 

0. 

0. 

M3 

9 

0. 

0. 

0. 

M3 

6 

0. 

0. 

0. 

M3 

7 

0. 

.100009*02 

.1000CE400 

M3 

a 

0. 

0. 

C. 

M3 

9 

0. 

0. 

0. 

M3 

10 

0. 

0. 

0. 

M3 

11 

.29900E-02 

-.9180CE*00 

.103606-02 

M3 

12 

0. 

0. 

0. 

M3 

13 

0. 

.10000E*01 

.2000CE-01 

M3 

14 

0. 

0. 

0. 

M3 

19 

0. 

0. 

.lOOOCE-02 

MS 

16 

0. 

0. 

0. 

M3 

11 

0. 

0. 

0. 
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L0«0-SET 


2 



SURFACE NODE REF. COORDINATES (INI RON 


NO. 

NO. 

BS/X 

BBL/y 

HL/Z 

N1 H2 

H3 

NO. 

1 

1 

900.0 

900.0 

2B0.0 


TX 

1. 

1 

1 

900.C 

900.0 

280.0 


Ty 

2. 

1 

1 

900.0 

900.0 

280.0 


TZ 

3. 

1 

1 

900.0 

900.0 

280.0 



4. 

1 

1 

900.0 

900.0 

280.0 

Rlr 


9. 

1 

i 

900.0 

900.0 

280.0 

Ri 


6. 


SURFACE 

NOOE 

REF. 

COORDINATES (INI 



ROM 

NO. 

NO. 

es/x 

BBL /V 

NL/Z 

Ml 

M2 M3 

NO. 

2 

1 

1290.0 

ICOO.O 

200.0 


TZ 

7. 


2 

1390.0 

ICOO.O 

200.0 


TZ 

8. 


3 

1200.0 

800.0 

2C0.0 


TZ 

9. 


4 

1390.0 

800.0 

2CO.O 


TZ 

10. 


9 

1190.0 

AOO.O 

20C.0 


TZ 

11. 


4 

1300.0 

600.0 

2C0.0 


TZ 

12. 


T 

1100.0 

490.0 

200.0 


TZ 

13. 


8 

1290.0 

490.0 

2C0.0 


TZ 

14. 


9 

1000.0 

300.0 

200.0 


TZ 

19. 


10 

1290.0 

300.0 

200.0 


TZ 

16. 


11 

900.0 

190.0 

2C0.0 


TZ 

17. 


12 

1200.0 

190.0 

2C0.0 


TZ 

18. 

'2 

1 

1290.0 

ICOO.O 

200.0 

TZ 


19. 


2 

1390.0 

ICOO.O 

200.0 

TZ 


20. 


3 

1200.0 

800.0 

200.0 

TZ 


21. 


4 

1390.0 

800.0 

200.0 

TZ 


22. 


9 

1190.0 

600.0 

2C0.0 

TZ 


23. 


6 

1300.0 

600.0 

200.0 

TZ 


24. 


T 

IIOO.O 

490.0 

200.0 

TZ 


29. 


8 

1290.0 

490.0 

2C0.0 

TZ 


26. 


9 

1000.0 

3C0.0 

200.0 

TZ 


27. 


10 

1290.0 

300.0 

200.0 

TZ 


28. 


112 


11 

900.0 

190.0 

2CO.O 

T2 

29. 

12 

1200.0 

190.0 

200.0 

T2 

30. 


SURFACE 

NODE 

REF. 

COORDINATES (INI 


RCh 

NO. 

NO. 

BS/X 

BSL/V 

Ml/2 

N1 N2 N3 

NO. 

3 

1 

190.0 

0.0 

290.0 

T2 

31. 


9 

1000.0 

0.0 

200.0 

T2 

32. 


SURFACE 

NOOE 

REF. 

COOROINATES (INI 



Rob 

NO. 

NO. 

BS/X 

BBL/V 

M./Z 

Nt M2 

M3 

NO. 

3 

1 

290.0 

0.0 

200.0 


TZ 

33. 


3 

400.0 

0.0 

200.0 


T2 

34. 


4 

S90.0 

0.0 

200.0 


TZ 

39. 


7 

IBOO.O 

0.0 

200.0 


TZ 

34. 



2090.0 

0.0 

210.0 


TZ 

37. 


9 

2290.0 

0.0 

230.0 


TZ 

38. 


10 

2490.0 

0.0 

290.0 


TZ 

39. 

3 

1 

290.0 

0.0 

200.0 

MV 


40. 


3 

400.0 

0.0 

200.0 

MV 


41. 


4 

890.0 

0.0 

200.0 

MV 


42. 


7 

1800.0 

0.0 

200.0 

MV 


43. 



2090.0 

0.0 

210.0 

MV 


44. 


9 

2290.0 

0.0 

230.0 

MV 


49. 


10 

2490.0 

0.0 

290.0 

MV 


44. 


SURFACE 

NOOE 

REF. 

COOROINATES (INI 




RCH 

ROTATION 

ANGLES (OEGI 

NO. 

NO. 

BS/X 

BBL/V 

ML/Z 

Ml 

HZ 

M3 

N0.THETAX THETAV 

THETAZ 

4 

1 

2490.0 

90.0 

400.0 



TZ 

47. 

0.0 

0.0 

”4.0 


3 

2470.0 

300.0 

400.0 



TZ 

4B. 

0.0 

0.0 

“4.0 

4 

1 

2490.0 

90.0 

4CO.O 


RX 


49. 

0.0 

0.0 

"4.0 


3 

2470.0 

300.0 

400.0 


RX 


90. 

0.0 

0.0 

•4.0 

4 

1 

2490.0 

90.0 

4CO.O 


R« 


91. 

0.0 

0.0 

“4.0 


3 

2470.0 

300.0 

400.0 


RV 


92. 

0.0 

0.0 

“4.0 

4 

1 

2490.0 

90.0 

400.0 

RX 



93. 

0.0 

0.0 

“4.0 


3 

2470.0 

300.0 

400.0 

RX 



94. 

0.0 

0.0 

“4.0 


Ml 

1 0. 

0. 

0. 

Ml 

2 0. 

0. 

0. 

HI 

3 0. 

0. 

0. 

Ml 

4 0. 

0. 

0. 

HI 

9 0. 

0. 

0. 



Ml 


0. 

0. 

C. 

HI 

7 

0. 

0. 

0. 

HI 

a 

0 , 

0. 

0. 

HI 

« 

0 . 

0. 

0. 

HI 

10 

0. 

0. 

e. 

HI 

11 

0. 

0. 

0. 

HI 

12 

0. 

0. 

0. 

HI 

13 

0. 

0. 

0. 

HI 

14 

0. 

0. 

0. 

HI 

IS 

0. 

0. 

0. 

HI 

16 

0. 

0. 

0. 

HI 

IT 

0. 

0. 

0. 

HI 

la 

0. 

0. 

0. 

HI 

19 

.10000E«01 

.SC00CE*02 

.lOOOCE^Ol 

HI 

20 

.10000E«01 

.15000E«03 

.BOOOCEtOO 

HI 

21 

.10000£*01 

0. 

.BOOOCE^OO 

HI 

22 

. 10000 E*01 

.15000E»03 

.7000CE«00 

HI 

23 

.10000E*01 

>.5COOOE*02 

.5000CE»00 

HI 

24 

.10000E*01 

.1C000E»03 

.4000CE«00 

HI 

29 

.10000E*01 

~.10000E*03 

.30000E«00 

HI 

26 

.10000E*01 

.9COOOE*02 

.20C0CE«00 

HI 

27 

.10000E*01 

-.2000CE*03 

.lOCOCE+00 

HI 

26 

.lOOOOE^Ol 

.9COOOE«C2 

.lOOOCE + 00 

HI 

29 

.10000E»01 

-.30000E«03 

.1000CE»00 

HI 

30 

.lOOOOEOl 

0. 

-.10000E»00 

HI 

31 

0. 

0. 

0. 

HI 

32 

0. 

0. 

0. 

HI 

33 

0. 

0. 

0. 

HI 

34 

0. 

0. 

0. 

HI 

3S 

0. 

0. 

0. 

HI 

36 

0. 

0. 

0. 

HI 

37 

0. 

0. 

0. 

HI 

36 

0. 

0. 

0. 

HI 

39 

0. 

0. 

0. 

HI 

40 

0. 

0. 

c. 

Ht 

41 

0. 

0. 

0. 

HI 

42 

0. 

0. 

c. 

HI 

43 

0. 

0. 

0. 

HI 

44 

0. 

0. 

c. 

HI 

45 

0. 

0. 

0. 

HI 

44 

0. 

0. 

c. 

HI 

47 

0. 

0. 

c. 

HI 

46 

0. 

0. 

0. 

HI 

49 

0. 

0. 

0. 

HI 

50 

0. 

0. 

c. 

HI 

91 

0. 

0. 

0. 

HI 

92 

0. 

0. 

0. 

HI 

53 

0. 

-.69796E-01 

-.99756E-0Z 

HI 

94 

0. 

.T3010E-04 

-.313Z2E-01 


H2 

1 0. 

0. 

0. 

H2 

2 0. 

0. 

0. 

H2 

3 0. 

0. 

c. 

H2 

4 0. 

0. 

.lOOOCE-01 
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HZ 

S 

0. 

.100006*01 

.200006-01 

HZ 

« 

0. 

0. 

.100006-02 

HZ 

7 

0. 

0. 

0. 

HZ 

• 

0. 

0. 

C. 

HZ 

« 

0. 

0. 

0. 

HZ 

10 

0. 

0. 

0. 

HZ 

11 

0. 

0. 

0. 

HZ 

IZ 

0. 

0. 

0. 

HZ 

13 

0. 

0. 

0. 

HZ 

14 

0. 

0. 

.0. 

HZ 

IS 

0. 

0. 

0. 

HZ 

1« 

0. 

0. 

0. 

HZ 

IT 

0. 

0^ 

0. 

HZ 

11 

0. 

0. 

0. 

HZ 

IS 

0. 

0. 

0. 

HZ 

30 

0. 

0. 

0. 

HZ 

31 

0. 

0. 

0. 

HZ 

ZZ 

0. 

0. 

0. 

HZ 

23 

0. 

0. 

0. 

HZ 

Z4 

0. 

0. 

0. 

HZ 

zs 

0. 

0. 

0. 

HZ 

z« 

0. 

0. 

0. 

HZ 

27 

0. 

0. 

0. 

HZ 

zo 

0. 

0. 

0. 

HZ 

zs 

0. 

0. 

0. 

HZ 

30 

0. 

0. 

0. 

HZ 

31 

0. 

0. 

0. 

HZ 

3Z 

0. 

0. 

0. 

HZ 

33 

0. 

0. 

0. 

HZ 

34 

0. 

0. 

0. 

HZ 

3S 

0. 

0. 

0. 

HZ 

36 

0. 

0. 

0. 

HZ 

37 

0. 

0. 

0. 

HZ 

30 

0. 

0. 

0. 

HZ 

3S 

0. 

0. 

0. 

HZ 

40 

0. 

.100006*01 

-.300006-01 

HZ 

41 

0. 

.100006*01 

-.100006-01 

HZ 

42 

0. 

.100006*01 

0. 

HZ 

43 

0. 

.100006*01 

.100006- 01 

HZ 

44 

0. 

.100006*01 

.200006-01 

HZ 

4S 

0. 

.100006*01 

.400006 -01 

HZ 

46 

0. 

.100006*01 

.700006-01 

HZ 

47 

0. 

0. 

0. 

HZ 

4S 

0. 

0. 

C. 

HZ 

4S 

0. 

>.647566-01 

-.99 7566 -02 

HZ 

SO 

0. 

.730106-04 

-.313226-01 

HZ 

SI 

0. 

.9S7S66*00 

-.697566 -03 

HZ 

S2 

0. 

.100246*01 

.176996-01 

HZ 

S3 

0. 

0. 

0. 

HZ 

S4 

0. 

0. 

0. 

H3 

1 

0. 

. 100006*02 

.100006*00 

H3 

2 

0. 

0. 

.120006*01 

H3 

3 

.10000E»01 

-.300006*03 

-.160006*01 


H3 
H3 
Hi 
Hi 
Hi 
Hi 
Hi 

Hi 

Hi 14 

Hi ,7 

"3 li 

Hi - 


* 0 . 

i *• 

« 0 . 


0 . 

0 . 

0 . 


io 

u 

ti 

ti 


w. 

•looaoc+oi ®* 
•tOOOOf*0| 

.10000E*0l ,*}^0O0e*03 
•WOOOf^Oi ‘ifooct^oj 

I? 'iS°°oe*03 

IM •*®000f»0l •*®0 ‘>OE»o? 

« « 0. 0. 

Hi in “• ‘ 

Jr 0- 0. c 

0 . 

0 . 


«3 

« - 
M3 
Hi 
Hi 
Hi 

Hi ,, 
30 


M3 

M3 


/4- ol- 

0. 
f* 0. 
0 . 
0 . 


0 . 

0 . 

0 . 

0. 

0 . 


3i 

33 


0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 

0 . 


:;Sns:si 

•jooooe*;^ 
7fJ®oe*oo 

*.a5®<**00 

**®foce»oo 

^•ioooce^oS 

c. 


M3 3/ 

Ma _ * 


Hi 

Hi 


34 

35 

30 

37 

30 

• 

|0 0 . 

0 . 

♦3 0. 
M 0. 

?* ®* 

*0 0 . 


0. 0, 

*10000f*oi ®* ?* 

•10000f*o| r*iJ?500f*04 

•WOOOC«01 .•if°00f»03 

• }f!°00t*Ot - 7J®®«E*03 

•10000f«0| .‘‘fOOflc^oj '?5®OCE»00 

• I0000f*03* • 330O0e*o/ . • f ®®OCE»00 

*fjf“®f*03 0**®“®*»00 

‘^oooilto} •,7®ooI*04 


0 . 

0 . 

0 . 

0 . 

0 . 

0 . 


•i0000f*0i 

♦ Woooe*oj ‘ii^°oe*o* 
„‘tfTooe*Si 


0 . 

0 , 

0 . 

0 . 

0 . 

0 . 


7- °®E*00 

o;Woooe»o® 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

•W00CE*00 

^.40000e»oo 


C. 

0 . 

0 . 

0 . 

C. 



LOAO-SET 


3 


SURFACE 

NODE 

REF. 

COOROINATES IINI 


ROM 

NO. 

NO. 

■SFX 

ML/V 

M./Z 

N1 M2 M3 

NO. 

2 

1 

1300.0 

ICCO.O 

200.0 

T2 

1. 


2 

1350. 0 

400.0 

2C0.0 

TZ 

2. 


» 

1000.0 

100.0 

200.0 

TZ 

3. 


4 

400.0 

130.0 

200.0 

rz 

4. 


S 

1000.0 

150.0 

200.0 

TZ 

5. 


N3 

1 -.100002*01 

-.10000E*03 

-. t04S4E*00 

NS 

2 -.10000E*01 

-.1S000E*03 

-. 7S31IE*00 

MS 

3 -.10000f*01 

.2000CE*03 

.13t5EE*0e 

M3 

4 -.100002*01 

.3flOOOE*03 

-.10000E*00 

M3 

5 -.100002*01 

.2C000E*03 

.•2320E-O1 


LOAI>>SET 


4 


PANEL LOADS I NORMAL 10 SLR FACE I 

LOAD-SET* 4 

FREQUENCIES • 0.000 2.000 


AERO 

SURFACE ROM NODE REF.COORDINA1ESI INI 


NO. 

NO. 

BS/X 8BL 

n hL/l 

1. 

1 

1200.0 

lOCO.O 

2C0.0 

2. 

2 

12SC.0 

1000.0 

200.0 

3. 

3 

1310.0 

ICOO.O 

200.0 

4. 

4 

1380.0 

lOCO.O 

200.0 

S. 

S 

IISO.O 

750.0 

200.0 


6 

1210.0 

750.0 

200.0 

T. 

7 

1280.0 

7S0.0 

200.0 

8. 

a 

1330.0 

750.0 

200.0 

9. 

9 

1080.0 

590.0 

200.0 

10. 

10 

IISO.O 

590.0 

200.0 

11. 

11 

1240.0 

590.0 

200.0 

12. 

12 

1300.0 

590.0 

200.0 

13. 

13 

1000.0 

400.0 

200.0 

14. 

14 

1100.0 

400.0 

200.0 

IS. 

IS 

1200.0 

4CQ.0 

200.0 

16. 

16 

1270.0 

400.0 

200.0 

IT. 

17 

900.0 

200.0 

200.0 

IS. 

18 

1030.0 

200.0 

200.0 

19. 

19 

1120.0 

200.0 

200.0 

20. 

20 

1200.0 

200.0 

200.0 

21. 

1 

24SC.0 

450.0 

400.0 

22. 

2 

2490.0 

450.0 

400.0 

23. 

3 

2930.0 

450.0 

400.0 

24. 

S 

2470.0 

150.0 

400.0 

2S. 

7 

2380.0 

40.0 

400.0 

26. 

8 

2430.0 

40.0 

400.0 

2T. 

9 

2500.0 

40.0 

400.0 


FREQUENCY 


l- 


0.000 CPS 


PLDS OPTION 


OAR MATRIX NERCEO 

TAPE tOAO 
ORDER SEO. 


1 

1. 

0. 

.1COC0E404 

-.10000E»02 

2 

2. 

0. 

.2CCC0E»C3 

-.70000E+01 

3 

3. 

0. 

.1Q000E*03 

-.70000E»01 

4 

4. 

0. 

.9COCCE402 

-.900C0E«01 

S 

9. 

0. 

.1COQ0E4O4 

-.80000E+01 

* 

6. 

0. 

.2C0C0E»03 

-.60000E+01 

7 

7. 

0. 

.lCOCOE+03 

-.40000E»01 

• 

B. 

0. 

.5COC0E«’O2 

-.20000E401 

• 

9. 

0. 

.aCOOOE+03 

-.aooocE^oi 

10 

10. 

0. 

.TOOCOE^-OS 

-.6C0C0E»C1 

11 

11. 

0. 

.TCOCCEt-OS 

-.40000E>C1 

12 

12. 

0. 

.1C000E403 

-.20000E401 

13 

13. 

0. 

.120C0E404 

-.TOOOOEtOl 

14 

14. 

0. 

• .6aOCOE«03 

-.40000E«’01 

15 

19. 

0. 

.30000E+03 

-.20000E»01 

14 

14. 

0. 

.1C000E403 

-.lOOOOEtCl 

IT 

IT. 

0. 

.14000E404 

0. 

U 

IS. 

0. 

.TCOCOE+03 

0. 

19 

19. 

0. 

.4COOOE*03 

0. 

20 

2J. 

0. 

,2cocce*a3 

0. 

21 

21. 

0. 

.2C0G0E403 

.100COE401 

22 

22. 

0. 

.9COCOE402 

.EOOOCE^CO 

23 

23. 

0. 

.1C0C0E*02 

.lOOOOEtOO 

24 

24. 

0. 

.lOOOOE+03 

.40000E4’00 

29 

29. 

0. 

.3COOOE403 

.20000E*C1 

26 

26. 

0. 

.100COE403 

.4aoooE+oo 

27 

27. 

0. 

.2COCOE402 

.20000E*00 


119 


L 



MS BAR MATRIX MERGED FREQUENCY 1- 0.000 CFS FLOS 

TAFE LOAD 
ORDER SEQ. 


1 

1. 

.lOOOOEtOl 

0. 

.100C0E»00 

2 

2. 

.2000CE«0C 

0. 

.TOOOOE-Ol 

J 

3. 

. taooo£*oa 

0. 

.700CCE-C1 

4 

4. 

.SOOOCE-01 

0. 

.SOOOOE-01 

S 

S. 

.10000E40I 

-.IC000E«00 

.aooooE-01 

* 

4. 

.20000E400 

-.1C0C0E«01 

.60000E-01 

7 

7. 

.10000E400 

0. 

.40000E-01 

a 

a. 

.S0O00E>01 

0. 

.200COE-01 

4 

4. 

.aooooEtoo 

-.20000E-01 

.aooooE-ci 

10 

10. 

.2OO0OE*OC 

-.ICOOOE-Ol 

.40000E-Ot 

11 

11. 

.20000E«00 

0. 

.40000E-ai 

12 

12. 

.IOOOOE»OC 

.IC0C0E*01 

.20000E-01 

13 

13. 

.I2000E*01 

-.4C0C0E-01 

.40000E-01 

14 

14. 

,40000E«Oa 

-.2COCOE>01 

.40000E-01 

IS 

IS. 

.30000E»OC 

-.2COO0E-O1 

.20000E-01 

14 

14. 

.lOOOOEtOO 

-.lCOCOE-01 

.lOOOOE-01 

17 

IT. 

.140CCE»01 

-.SCOCCE-01 

0. 

IS 

IS. 

.70000E*00 

-.3C0CCE-01 

0. 

19 

14. 

.40000E«00 

>.2CCC0E-0l 

0. 

20 

20. 

.2OOOCE«0C 

-.1CCOCE>01 

0. 

21 

21. 

.2000CE4-00 

.2CO00E*O3 

-.SOOOOE»C2 

22 

22. 

.SOOOCE-Ol 

.5CO0OE*O2 

-.10000E«^02 

23 

23. 

.lOOCCE-01 

.1C0C0E»02 

-.S0000E*01 

24 

24. 

.lOOOOEtOO 

.IOOOOE403 

-.SOOOOE*01 

2S 

2S. 

.30030E4’OC 

.300C0E403 

-.20000E4’02 

24 

24. 

.10000E«OC 

.100C0E«03 

0. 

27 

27. 

.20000E-01 

.20000E«’02 

.SOOOOEtCl 


OFTION 


120 



C3 BAR lOR PHI-riLOA BARI NATRtR MERGED PREGUENCY !■ 0.000 CPS PLOS OPriON 

YAPE LOAD 
ORDER SEQ. 



.1000E*01 0. 

,zoooe*oo 0 . 
.tOODE«O0 0. 
.SOOOE-01 C. 
.lOOOEtOl 0. 
.2000E«00 0. 
.1000E«00 0. 
.9000E-01 0. 
.B00DE«0D 0. 
.200DEA00 0. 
.2000E»00 0. 
.1000E»00 0. 
.1200E«01 0. 
.A00DE*O0 0. 
.300CE»00 0. 
.1000E*00 0. 
.1400E*01 0. 
.TOOOEAOO 0. 
.A000E«00 0. 
.2000E*00 0. 
.2000E«00 0. 
.SOOOE-01 0. 
.lOOOE-01 0. 
.1000E«00 0. 
.3000E*0a 0. 
.1000E400 0. 
.2000E-01 0. 
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H4 tAR MATRIX MERGED FREOUENCV 2- 2.000 CPS PLOS 

TAPE LOAD 
ORDER SEO. 


1 

1. 

0. 

.1C0C0E403 

-.10000E40I 

2 

2. 

0. 

.2CCC0E402 

-.70000E*00 

3 

3. 

0. 

.100COE402 

-.70000E400 

A 

4. 

0. 

.500CCE401 

-.500CCE400 

9 

5. 

0. 

.lOOCOE^OS 

-.saoccE*oo 

A 

4. 

0. 

.20000E402 

-.400COE400 

T 

7. 

0. 

.1CCC0E402 

-.40000E400 

R 

• 

0. 

.5C0C0E*01 

-.20000E*00 

* 

9. 

0. 

.RCOCOE402 

-.SOO0CE40O 

10 

10. 

0. 

.200COE402 

-.40000E400 

11 

11. 

0. 

.20000E402 

-.400COE*00 

12 

12. 

0. 

.1C0C0E«02 

•.20000E400 

13 

13. 

0. 

.120C0E403 

-.400Q0E»00 

14 

14. 

0. 

.600C0E+O2 

-.400C0EP00 

15 

15. 

0. 

.3C0C0E402 

-.20000E400 

1* 

14. 

0. 

.10CCOE402 

-.10000E400 

17 

17. 

0. 

.140COE403 

0. 

IR 

IS. 

0. 

.7COOOE+02 

0. 

19 

19. 

0. 

.4C0C0E402 

0. 

20 

20. 

0. 

.2C0C0E«02 

0. 

21 

21. 

0. 

.2C0COE*03 

.10000E»01 

22 

22. 

0. 

.500COE402 

.2000CEtOO 

23 

23. 

0. 

.1COOOE402 

.10000E«00 

24 

24. 

0. 

.1C0C0E*03 

.40000E*00 

25 

25. 

0. 

.30000E403 

.10000E4 01 

24 

24. 

0. 

.1CO0OE4O3 

.20000E400 

2T 

27. 

0. 

.20000E402 

.10000E400 


OPTION 



MEOUENCV 2> 2.000 CPS 


PLOS OPTION 


I 


MS BAR MATRIX MERCED 
TAPE LOAD 
CRDER SEO. 


1 

1. 

.1000CE«01 

0. 

.lOOOOE^OO 

2 

2. 

.2000CEtOC 

0. 

.700COE-01 

3 

3. 

.10000E400 

0. 

.700C0E-01 

4 

4. 

.SOOOCE-01 

0. 

.SOOCOE-Ol 

S 

5. 

.10000E401 

-.ICOOOEtOO 

.800COE-01 

6 

6. 

.20000E»Oa 

-.lCOOOE-01 

.60000E-C1 

7 

7. 

.1000CE40C 

0. 

.40000E-01 

a 

8. 

.SOOOCE-01 

0. 

.20000E-01 

9 

9. 

.SOO0CE4O0 

-.200COE-01 

.80000E-01 

1C 

10. 

.2000CE«0C 

-.lOOCOE-01 

.60000E-01 

11 

11. 

.20000E«00 

0. 

.4Q0C0E*01 

12 

12. 

.10000E400 

.lOOOCE-01 

.20000E-01 

13 

13. 

.12000E601 

-.4COCOE-01 

.40000E-C1 

14 

14. 

.60000E»0e 

-.2C0C0E-01 

.40000E-01 

15 

15. 

.30000E+00 

-.2CCC0E-C1 

.20000E-01 

16 

16. 

.10000E600 

-.lCOflOE-01 

.IOOOOE-01 

17 

17. 

.14000E401 

-.SCOOOE-01 

0. 

18 

18. 

.70000E*00 

-.3C0C0E-01 

0. 

19 

19. 

.40000E«0C 

-.2C0C0E-01 

0. 

20 

20. 

.20000E«00 

-.ICOCOE-Ol 

0. 

21 

21. 

.200008 *C0 

.1C0C0E403 

-.SOOCOE«C1 

22 

22. 

.SOOOOE-01 

.2C0C0E402 

-.10000E401 

23 

23. 

.lOOOOE-01 

.10000E«02 

-.SO000E60O 

24 

24. 

.1000CE*00 

.20000E*02 

-.5OOC0E4CO 

25 

25. 

.30000E*oe 

.1COOOE403 

-.20000E»01 

26 

26. 

.10000E40C 

.2CCC0E402 

0. 

27 

27. 

.20000E>01 

.1COOOE«02 

.SOOOOE*00 
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FREOUENCr !■ 2.000 CPS 


PL OS OPTION 


C3 BAR (OR PHI-TILOR BAR) MATRIX MERGED 


tape 

LOAD 



ORDER 

SEO. 



1 

1. 

.1000E401 

.10C0E*00 

2 

2. 

.20006400 

.IOCCE*CO 

S 

3. 

.1000E400 

.10COE*CO 

4 

4. 

.50006-01 

.10006*00 

5 

S. 

.10006*01 

.20006*00 

6 

6. 

.20006*00 

.20006*00 

7 

7. 

.10006*00 

.20006*00 

B 

B. 

.90006-01 

.20006*00 

9 

9. 

.S000E*00 

. 10006*00 

10 

10. 

.20006*00 

.10006*00 

11 

11. 

.20006*00 

. 10006*00 

12 

12. 

. 10006*00 

.10006*00 

13 

13. 

.12006*01 

0. 

14 

14. 

.60006*00 

c. 

IS 

15. 

.30006*00 

0. 

16 


.1000E*00 

0. 

17 

17. 

.14006*01 

0. 

IS 

IS. 

.70006*00 

0. 

19 

19. 

.400C6*00 

0. 

20 

20. 

.20006*00 

0. 

21 

21. 

.20006*00 

-.20006*00 

22 

22. 

.50006-01 

-.50006-01 

23 

23. 

.IOOOE-01 

-.10006-01 

24 

24. 

.10006*00 

-.10006*00 

25 

25. 

.30006*00 

-.30006*00 

26 

26. 

.10006*00 

-.10006*00 

27 

27. 

.2000E-01 

-.20006-01 
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LQAD'SE? 


9 


ET PANEL LOADS I NORMAL TO SURFACE I 
LOAD-SETo 

FREQJENCIES °> OoOOO 2.0C0 


S 


STRUCTo 

SURFACE ROW NODE REF COORDINATESUNi 
NU< NOa NOo BS/8 BBL/r NL/Z 



1 

129CcO 

1000.0 

200.0 

3o 

2 

1330c0 

1000.0 

200.0 

3. 

3 

IZOCcO 

SOO.O 

200.0 


A 

135CcO 

8C0.0 

200.0 

9. 

9 

llSOcO 

6CO.O 

200.0 

&o 

6 

ISOOoO 

&C0.0 

200.0 

Te 

T 

llOCcO 

490.0 

2C0.0 

Oo 

0 

1290c0 

490.0 

200.0 

9c 

9 

1000.0 

300.0 

200.0 

tOc 

10 

129C.0 

300.0 

200.0 

lie 

11 

9CO.O 

190.0 

200.0 

I2c 

12 

1200.0 

190.0 

200.0 

I3c 

1 

2A3Ce0 

30.0 

400.0 

I9e 

3 

246 CcO 

190.0 

400.0 

l»c 

3 

24TO.O 

300.0 

400.0 


N9 BAA MERGED HATHIlt 

LOAD SEQo 


1 

0. 

0. 

0. 

2 

0. 

c. 

0. 

3 

0. 

c. 

0. 

4 

0. 

c. 

0. 

3 

0. 

c. 

0. 

6 

0. 

c. 

0. 

1 

0. 

0.- 

0. 

0 

0. 

c. 

0. 

9 

0. 

c. 

0. 

10 

0. 

c. 

0. 

11 

Oo 

0. 

0. 

12 

0. 

c. 

0. 

13 

.lOOOOE^Ol 

.129C0E«04 

.TOOOOE400 

14 

.2COOOE»01 

.25200E«04 

.12000E4C1 

19 

.30000E»01 

.3aiC0E«04 

.12000E4C1 
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F14 BA0 MATfttX HERCEO FRECUENCV 1 « 0«000 CPS NPLDS OPTION 

TAPE LOAD 
QflDEP SE4. 


1 

1. 

0. 

,6A000E»01 

-.23iaCE*C0 

2 

2. 

0. 

.25945£*C1 

-.19994E*00 

3 

3. 

0. 

,12S26E»02 

-.19146E4'00 

A 

A. 

0. 

.15Q35E»C1 

-.70303E-01 

5 

5. 

0. 

.T3734E»01 

•.14552EP00 

& 

6. 

0. 

.27012E«01 

-.54163E-01 

7 

7. 

0. 

,1166DE*0Z 

-.U277£<-00 

0 

S. 

0. 

.34522E«ai 

-.39578E-C1 

9 

9. 

0. 

.lS27aE*a2 

-.12771E-0I 

lO 

10., 

0. 

.I7227Ea<U 

-.47207E-02 

11 

11. 

a. 

.l$9a8E*02 

.17610E-0I 

12 

12. 

0. 

.217ClE»ai 

.12299E-01 

13 

13. 

0. 

.630C0Et03 

.390COE*01 

lA 

14. 

0. 

.34QOaE*03 

.19SOOEAC1 

1$ 

15. 

0. 

ol30a0EA03 

.65000E«00 
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H5 Bfip nATKIX OERCEO FBEflUENCV 1= 0„000 CPS NPLDS OPTION 

TAPE LOAD 
OHDEQ SEO. 


1 

le 

efi6000E-02 

-.066S1E>1T 

.23100E-02 

2 

2. 

.3594S£*02 

.41940E«a5 

.199S7E-02 

3 

3. 

.I2S2AE-01 

~.l23C7E-02 

.191376-02 

4 

4 . 

olS075E~02 

-,53262E'C4 

o70237E-C3 

5 

5. 

.J3734C-02 

-.4S726E-03 

.163936-02 

6 

6. 

.27072E-02 

.1924$.E-03 

.361286-03 

7 

7 , 

.llliEOE-Ol 

-.31537E-03 

.11364E-02 

0 

0. 

.34S22E->02 

->.U34aE>03 

.39573E-03 

9 

9o 

.ia270E^Cl 

->.6113AE-03 

.35871E-03 

10 

10. 

.17227E-02 

-.1S914E-03 

.48509E-04 

11 

lU 

.I4988E-0I 

-.765176-03 

-.219236-03 

12 

12. 

.2l7aiE-02 

-.13626E-01 

-.1317BE-C3 

13 

11 . 

.61000E«'OC 

.03CC0E»03 

-.32500EP02 

14 

14. 

.34000E»0C 

.34000EP03 

-.50000EP02 

15 

15. 

.iaoooE«oo 

.13000EP03 

-.32500EP02 
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FREOUENCr 


1 > 


0.300 CPS 


NPLDS DP non 


CO Bill (DP Phl-TILDA BAP > 
TAPE LOAO 
ORO£A SEO. 


HATRIX MERCED 


1 

1. 

.«600E-OZ 

0. 

z 

Z. 

.Z594E-02 

0. 

3 

3. 

.12S3E-01 

0. 

4> 

A. 

. 1507E-0Z 

0. 

5 

5. 

.7373E-0Z 

0. 

6 

b. 

.Z707E-0Z 

c. 

7 

7. 

. 1166E-01 

0. 

8 

8. 

.36SZE-0Z 

c. 

9 

9. 

.18Z7E-01 

0. 

10 

10. 

.17Z3E-0Z 

0. 

u 

11. 

.1999E-01 

0. 

IZ 

IZ. 

-217OE-0Z 

0. 

13 

13. 

.«300E*00 

0. 

lA 

lA. 

.BAOOE^OO 

0. 

13 

IS. 

.1300EA00 

0. 


128 



R(M 

LOAD-SET 6 

SURFACE LO&O REF CQaftDInaTES 

(INI 

aais oaiErtTATion 

IDEGI LOCO 


KOb 

HOb 

SEQb 

05/iI 

6BL/V 

13L/2 

THETC 

THET7 

THETZ 

TVPE 

1 

1 

I^Ob 

&o 

llOCoO 

SCOoO 

SCOoO 

OoO 

OoO 

0.0 

v» 

2 


lo 
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ROU 

NO. 


M3 e«R MERCED MATRIX COEFEICIENTS NUMBER OF MODES* 

1 0. .23900Et02 .25900£*CO 

2 0. C. ,25900E*C0 

3 .2S900EAOL ~.S18C0E'>C3 .10360E*01 

4 0. C. -.25900E»C1 

3 -.2S900E*03 .52836E*0S -.8IOAOE«C2 

6 0. 0. .26936E4 02 

T -.087a6E-Ol -.S983NE*C1 -.S67A6E-01 

8 0. .2S9CCE<'a2 -.259008*00 

9 -.2A39NE»00 --.15591E + 00 

10 0. C. .25900£i-0a 

11 .29£72£*01. .199968*03 .1<969E«C1 

12 .128Nb£*02 -.81585E*03 .3e047E*a 

13 .762026*03 .668A5E*05 .56153E*C3 

16 .38526E*06 ~.97669E*C3 .179036*06 

15 -.66560£*02 .TSBSTE'AQA -.666066*02 

16 .111116*06 .162266*06 .666466*03 

17 .176966*02 .262626*04 .113436*02 

la 0. -.906506*04 -.76666E*02 

19 0. C. 0. 

20 0. C. 0. 

21 .692106*00 .612:66*03 .279726*00 

22 .906506*02 .116606*06 .351726*02 

23 .725206*01 .926896*06 .295266*01 

26 0. 0. 0. 

25 0. 0. -.621606*00 

26 0. C. 0. 

27 .220156*01 -.16698£*C6 .375556*00 

28 .525776*01 .505836*06 .156956*01 

29 .113666*02 .641806*06 .11810E*0i 

30 0. 0. 0. 

31 -.589236*03 .420126*06 -.131836*03 

32 .199636*06 .224396*07 .106266*06 

33 .659936*06 .579696*07 .193356*06 

36 0. 0. 0. 


3 


130 



Rcm 
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HS BAR HATRIX NEfiCED FREQUENCY lo Q .000 CFS 
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RCtd 
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us BitR (4RTOIS UetCEO 
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0. 
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0. 
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0. 
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0. 

21 
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olOOOOEoOl 

20 

.22060E401 

«70013Efr03 
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2*3 

<.22020E40L 

o79ClSE4C3 

-<.22000E»02 

30 

0. 

0. 

Oo 

31 

-.160Q0EO03 

.84000E»C3 

~o4a000E»03 

32 

.13S39E404 
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-oS4S70E405 

33 

o27063E*04 

eO6074E»C6 

-o20310£40S 

34 

0. 

0. 

0. 


C} (OR PHI'TllDA OAR) MATRIA PfRCED fRECUENCV 2« 2.000 CPS 

ROW 

NO. 


i 

0. 

0. 

2 

0. 

0. 
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.2000E'»C2 

0. 

A 

0. 

0. 

S 
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0. 
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0. 

0. 

T 
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6 

0. 

0. 
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0. 
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22 
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~.ZBOOBtOl 

.280CE»C: 

24 

0. 

0. 

25 

0. 

0. 

26 

0. 

0. 

27 

.5000E400 

.30CCE*00 

26 

.2206E401 

-.2U4E*0l 

29 

.220ZE401 

•'.219BE4C1 

30 

0. 

0. 

31 

-.1600E403 

--160CE+03 

32 

.13S4E*C4 

-.1377£»« 

33 

.270&E404 

-.269SE«04 

34 

0. 

0. 
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